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PREFACE 


This  report  presents  the  research  accomplishments  for  FY87  (1  October  1986  to  30 
September  1987)  of  the  research  investigation  entitled  Theoretical  Investigation  of  Three- 
Dimensional  Shock-Wave  Turbulent  Boundary  Layer  Interactions". 

The  research  has  benefited  from  the  assistance  of  several  individuals,  including  Drs. 
James  Wilson  and  Jim  McMichael  (Air  Force  Office  of  Scientific  Research),  Dr.  C.  Horstman 
(NASA  Ames  Research  Center),  and  Dr.  James  Keller  and  Mr.  Manuel  Salas  (NASA  Langley 
Research  Center).  The  interactions  with  S.  Bogdonoff,  D.  Dolling,  A.  Ketchum,  and  G.  Settles 
are  acknowledged. 
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Section  I.  Objectives 


The  principal  objectives  of  the  research  program  are  : 

1.  Develop  and  Validate  Theoretical  Model(s)  for  3-D  Shock  Wave  •  Turbulent  Boundary 
Layer  Interaction 

The  understanding  of  complex  3-D  shock  wave-turbulent  boundary  layer  interactions  ("3- 
D  turbulent  interactions")  requires  the  development  and  validation  of  accurate  theoretical 
models.  A  wide  range  of  theoretical  approaches  have  been  employed,  extending  from 
simplified  control  volume  analyses  (Paynter  1980)*  to  the  Reynolds-averaged  three- 
dimensional  compressible  Navier-Stokes  equations  (Horstman  and  Hung  1979;  Knight 
1984a,  1984b,  1985a;  Knight  et  al  1986)  with  turbulence  incorporated  through  a  tur¬ 
bulent  eddy  viscosity  model.  The  present  research  effort  has  utilized  the  algebraic  tur¬ 
bulent  eddy  viscosity  of  Baldwin  and  Lomax  (1978).  The  equations  are  solved  numeri¬ 
cally  using  a  hybrid  explicit-implicit  numerical  algorithm  developed  by  the  present  investi¬ 
gator  (Knight  1984a).  The  validation  is  achieved  by  comparison  of  computed  and  experi¬ 
mental  results  for  3-D  turbulent  interactions  in  simplified  geometries.  This  validation  has 
been  a  continuous  element  of  the  research. 


2.  Determine  the  Physical  Structure  of  3-D  Turbulent  Interactions  for  Selected 
Geometries 

A  second  objective  is  the  determination  of  the  physical  flowfield  structure  of  3-D  tur¬ 
bulent  interactions  in  simplified  geometries  (e.g.,  3-D  sharp  fin  [Fig.  1]  and  3-D  swept 
compression  corner  [Fig.  2]).  This  understanding  is  naturally  linked  to  the  development 
and  verification  of  accurate  theoretical  model(s)  for  3-D  turbulent  interactions,  and  as 
such  requires  a  coordinated  effort  between  theory  and  experiment.  During  FY85**  the 
collaboration  of  C.  C.  Horstman  (NASA  Ames  Research),  B.  Shapey  and  S.  Bogdonoff 
(Princeton)  and  the  author  led  to  an  understanding  of  the  flowfield  structure  of  the  3-D 
sharp  fin  configuration  (Knight  et  al  1987a).  During  FY86,  a  similar  collaboration  led  to 
the  understanding  of  the  flowfield  structure  for  the  3-D  swept  compression  comer  at 
Mach  3  for  (a.  A.)  =  (24  ,  60)  deg.  During  FY87,  the  flowfield  structure  for  the  3-D 
swept  compression  comer  at  Mach  3  for  (a,  A.)  =  (24,  40)  deg  was  determined,  including 
a  detailed  quantitative  description  of  the  surfaces  of  separation  and  attachment. 


•  References  are  included  in  Section  V. 

**  Fiscal  Years  (FY)  refer  to  the  Federal  fiscal  year  (i.e.,  October  to  September). 
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3.  Investigate  Methods  for  Control  and  Modification  of  3-D  Turbulent  Interaction 
Flowfields 

The  understanding  of  the  flowfield  structure  of  3-D  turbulent  interactions  provides  the 
opportunity  for  investigation  of  methods  for  control  and  modification  of  the  interaction 
flowfields.  During  FY87  a  significant  effort  was  focused  on  determining  the  effect  of  sur¬ 
face  bleed  on  the  flowfield  structure  of  the  3-D  sharp  fin  interaction. 


4.  Develop  a  Unified  Understanding  of  the  Flowfield  Structure  of  3-D  Turbulent  Interac¬ 
tions 

A  fourth  objective  is  the  understanding  of  the  flowfield  structure  within  broad  families 
of  3-D  turbulent  interactions.  The  research  effort  to  date  has  focused  on  the  family  of 
dimensionless  3-D  shock  generators  (e.g.,  3-D  sharp  fin  and  3-D  swept  compression 
corner).  The  results  have  shown  that  the  flowfield  structure  is  similar  for  both 
configurations,  namely,  a  large  vortical  structure.  Further  research  is  in  progress  to  deter¬ 
mine  the  universality  of  this  structure  to  other  configurations. 


The  research  program  has  achieved  significant  progress  towards  these  objectives.  The 
present  report  details  the  progress  during  FY87  (1  October  1986  -  30  September  1987),  and 
describes  the  research  program  for  FY88  (1  October  1987  -  30  September  1988). 


Section  II.  Research  Accomplishments  for  FY87 

and 

Research  Program  for  FY88 


A.  Research  Accomplishments  for  FY87  (1  October  1986  -  30  September  1987) 

The  research  program  has  followed  the  basic  tasks  outlined  in  the  original  proposal  sub¬ 
mitted  to  Dr.  James  Wilson  on  4  April  1986.  Four  specific  research  tasks  were  outlined  for 
FY87,  namely,  i)  the  development  and  validation  of  theoretical  models  for  3-D  shock  wave- 
turbulent  boundary  layer  interaction,  ii)  the  determination  of  the  fiowfield  structure  of  3-D  tur¬ 
bulent  interactions  for  selected  geometries,  iii)  the  investigation  of  the  control  and  modification 
of  3-D  turbulent  interactions,  and  iv)  the  development  of  a  unified  understanding  of  the 
fiowfield  structure  of  3-D  turbulent  interactions.  The  progress  during  FY87  in  these  four  areas 
is  summarized  in  this  section.  Complete  details  of  the  research  activity  are  provided  in  the 
papers  presented  at  the  AIAA  26th  Aerospace  Sciences  Meeting  (see  Section  VI). 

The  overall  research  program  is  closely  coordinated  with  the  computational  research  of 
Dr.  C.  C.  Horstman  (NASA  Ames  Research  Center)  and  the  experimental  research  of  Prof.  S. 
Bogdonoff  (Princeton  University).  Drs.  Knight,  Horstman  and  Prof.  Bogdonoff  meet  fre¬ 
quently  to  select  specific  3-D  shock  wave-turbulent  boundary  layer  interactions  for  detailed 
study.  Separate,  independent  computations  are  then  performed  for  the  configuration  by  Drs. 
Knight  and  Horstman,  and  the  calculated  flowfields  are  compared  with  the  experimental  meas¬ 
urements  of  Bogdonoff.  This  close  collaboration  between  theory  and  experiment  is  a  critical 
element  of  the  success  of  the  program. 


Research  Task  No.  1  :  Develop  and  Validate  Theoretical  Models 

for  3-D  Shock  Wave-Turbulent  Boundary  Layer  Interactions 

The  theoretical  model  is  the  3-D  Reynolds-averaged  compressible  Navier-Stokes  equa¬ 
tions,  combined  with  a  model  for  the  Reynolds  stresses  and  turbulent  heat  flux.  The  Reynolds 
stresses  are  modelled  using  a  turbulent  eddy  viscosity;  the  turbulent  heat  flux  is  obtained  from 
the  turbulent  Reynolds  stresses  assuming  a  constant  turbulent  Prandtl  number  equal  to  0.9. 
Three  separate  turbulent  eddy  viscosity  models  are  employed,  namely,  i)  the  algebraic  model  of 
Baldwin  and  Lomax,  ii)  the  algebraic  model  of  Cebeci  and  Smith,  and  iii)  the  k-e  model  of 
Jones  and  Launder.  The  computations  are  performed  by  the  principal  investigator  (using  the 
Baldwin -Lorn ax  model)  and  Dr.  C.  C.  Horstman  (using  the  Cebeci-Smith  and  k-e  models) 
employing  the  CYBER  205  and  CRAY-X/MP,  respectively. 

The  theoretical  models  have  been  successfully  applied  to  the  prediction  of  the  3-D 
sharp  fin  interaction  (Knight  el  al  1987a).  During  FY86  the  model  was  successfully  applied  to 
the  3-D  swept  compression  corner  for  (a,  \)  =  (24,  60)  deg  at  Reg^  =  1.4xl05  and  9xl0') 
(Knight  et  al  1987b).  The  computed  flowfields  using  the  two  turbulence  models  were  compared 
with  the  experimental  data  from  the  Princeton  Gas  Dynamics  Laboratory  for  surface  pressure, 
and  boundary  layer  profiles  of  pitot  pressure  and  yaw  angle.  For  the  lower  Reynolds  number 
cases,  the  computed  surface  pressure  displayed  significant  disagreement  with  experiment. 
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although  specific  features  (e.g.,  upstream  influence  location,  peak  corner  pressure)  were  accu¬ 
rately  predicted.  The  reason  for  the  disagreement  is  not  clear,  and  is  the  subject  of  current 
study.  For  the  high  Reynolds  number  case,  the  computed  flowfields  showed  good  agreement 
with  experimental  measurements  of  surface  pressure,  and  boundary  layer  profiles  of  pitot  pres¬ 
sure  and  yaw  angle  (Knight  et  al  1987b). 

During  FY87  the  theoretical  models  were  applied  to  the  3-D  swept  compression  corner 
for  a  different  sweep  angle,  specifically,  the  (a,  X)  =  (24,  40)  deg  configuration,  at  Reg^  = 
2.6xl05  and  8.1x10  .  Computations  were  performed  by  Knight  using  the  Baldwin-Lomax 
model.  Calculations  were  performed  earlier  by  Horstman  using  the  Jones-Launder  (k-e)  and 
Cebeci-Smith  turbulence  models.  The  computed  flowfields  were  compared  with  the  experimen¬ 
tal  data  of  Settles  and  McKenzie  (Settles  et  al  1986)  and  Ketchum  and  Bogdonoff  (Knight  et  al 
1988).  The  experimental  data  include  surface  pressure  and  boundary  layer  profiles  of  pitot 
pressure  and  yaw  angle.  A  detailed  comparison  of  theory  and  experiment  is  presented  in 
Knight  et  al  (1988)*.  The  results  of  the  comparison  may  be  summarized  as  follows  : 

1.  Surface  Pressure 

The  calculated  surface  pressure  obtained  using  the  Baldwin-Lomax  model  underpredicts 
the  position  of  upstream  influence  (i.e.,  the  location  of  the  initial  pressure  rise).  The 
predictions  using  the  Cebeci-Smith  and  Jones-Launder  models  are  in  closer  agreement 
with  experiment;  however,  additional  grid  refinement  studies  are  needed  to  determine  the 
effect  of  grid  resolution  for  the  computations  using  the  Cebeci-Smith  and  Jones-Launder 
models.  The  Baldwin-Lomax  model  provides  an  accurate  prediction  of  the  plateau  pres¬ 
sure  in  the  vicinity  of  the  corner  line,  while  the  Cebeci-Smith  and  Jones-Launder  models 
overpredict  the  plateau  pressure. 

2.  Pitot  Pressure  and  Yaw  Angle 

The  computed  pitot  pressure  and  yaw  angle  profiles  are  observed  to  be  remarkably  insen¬ 
sitive  to  the  turbulence  model,  and  in  good  agreement  with  experimental  measurements. 
The  computed  yaw  angles  in  the  immediate  vicinity  of  the  surface  (i.e.,  the  lower  10%  of 
the  boundary  layer)  are  observed  to  be  sensitive  to  the  turbulence  model,  and  differences 
of  10%  -  15%  are  observed. 

3.  Surface  Visualization 

The  computed  surface  skin  friction  lines  exhibit  qualitative  agreement  with  the  experi¬ 
mental  kerosene-lampblack  surface  visualization.  The  calculated  line  of  coalescence 
occurs  downstream  of  the  experimental  position,  with  a  typical  angular  difference  of 
approximately  10%.  The  line  of  divergence  is  evident  on  the  compression  ramp  surface  in 
both  the  computation  and  experiment. 


*  Included  in  Section  VI. 


Research  Task  No.  2:  Determination  of  Physical  Structure 

of  3*D  Turbulent  Interactions  for  Specific  Geometries 

The  research  program  has  focused  on  the  development  of  flowfield  models  for  3-D  tur¬ 
bulent  interactions.  This  effort  is  divided  into  two  major  categories,  namely,  l)  the  develop¬ 
ment  of  sophisticated  flowfield  analysis  and  visualization  tools,  and  2)  the  determination  of 
flowfield  models  using  graphical  techniques.  These  efforts  are  described  below. 

Development  of  Sophisticated  Flowfield  Analysis  and  Visualization  Tools 

During  FY86  a  major  new  facility  was  established  in  the  College  of  Engineering  to  facil¬ 
itate  large  scale  scientific  computation.  This  facility  includes  a  Convex  C-l  minisupercomputer, 
two  Sun  3/180  fileservers  with  l.S  GByte  disk  storage,  nine  Sun  graphics  workstations  (includ¬ 
ing  two  color  workstations),  and  a  variety  of  hardcopy  devices  (color  thermal  hardcopy  and 
laser  printers  (4).  This  facility  was  developed  with  funds  received  from  i)  the  State  of  New 
Jersey’s  Commission  on  Science  and  Technology,  and  ii)  Rutgers  University.  Additional  sup¬ 
port  has  been  received  from  the  Air  Force  Office  of  Scientific  Research.  The  current  Director 
of  the  facility  is  Prof.  Abdel  Zebib,  Department  of  Mechanical  and  Aerospace  Engineering,  who 
succeeded  Prof.  Doyle  Knight  on  1  October  1987. 

A  significant  effort  has  focused  on  development  of  software  tools  for  the  visualization 
of  streamlines  in  order  to  ascertain  the  flowfield  structure.  The  following  tasks  have  been 
accomplished  during  FY87  : 

1.  Development  of  a  general  particle  tracing  code 

A  generalized  three  dimensional  particle  tracing  computer  code  was  completed  FY87.  The 
program  permits  a  variety  of  user-specified  interpolation  schemes  (linear  and  quadratic)  to 
integrate  the  motion  of  a  fluid  particle  within  an  arbitrary  grid  system.  A  complete 
description  of  the  algorithm  is  provided  in  Raufer  (1987).  The  program  was  successfully 
tested  against  a  previous  code,  modelled  after  the  particle  tracing  algorithm  of  Pieter  Bun- 
ing,  NASA  Ames  Research  Center  (Buning  1984).  The  code  has  been  extensively  util¬ 
ized  in  determining  the  flowfield  structure  of  the  3-D  swept  compression  corner  (see 
below). 

2.  Development  of  graphics  software  for  visualization  of  particle  traces 

A  sophisticated,  menu-driven  graphics  software  program  for  the  Sun  workstations  was 
completed  in  FY87.  The  program,  written  by  Dr.  Sandy  Walther  (Consultant),  provides  a 
convenient  interface  for  examining  particle  pathlines.  The  tiled  menus  permit  rapid  selec¬ 
tion  of  user-specified  visualization  parameters  (e.g.,  location  of  the  eye  of  the  observer, 
picture  plane,  etc.). 

3.  Development  of  Analysis  Code  for  Determining  Viscous  and  Inviscid  Effects 

A  research  effort  was  initiated  in  FY87  to  develop  a  computer  code  to  examine  the  contri¬ 
butions  to  the  rate  of  change  of  the  mean  kinetic  energy  ("mechanical  energy")  along  a 
flow  streamline.  The  mechanical  energy  equation  is 


D^Ujup/Dt  =  -(u-/p)5p/3xj  +  (ui/p)3xiy,3xj 

where  ViUjUj  is  the  mean  kinetic  energy  per  unit  mass,  D/Dt  =  3/3t  +  u-3/3x-, 
(Uj.U2.U2)  is  the  cartesian  velocity  vector,  (xj.x^.x^)  are  the  cartesian  coordinates,  p  Is 
the1  static  pressure,  p  is  the  density,  tjj  is  the  total  (viscous  plus  turbulent)  stress  tensor, 
and  the  Einstein  summation  convention  is  employed.  There  are  two  contributions  to  the 
rate  of  change  of  the  mean  kinetic  energy  per  unit  mass.  The  first  term,  -(uj/p)3p/3xj,  is 
an  in  viscid  effect  associated  with  compression  or  dilation.  The  second  term, 
(u,/p)3tij/3xj ,  represents  the  effect  of  the  viscous  and  turbulent  stresses.  A  series  of 
flow  streamlines  are  generated  by  integrating  the  motion  of  a  fluid  particle,  released  at  a 
selected  upstream  position,  through  the  computed  flowfield.  The  mean  kinetic  energy 
equation  is  integrated  simultaneously  along  the  streamline,  i.e., 

!4ujUj  =  ViuiUil0  +  /  {-(Uj/p)3p/3xj  +  (ui/p)3Ty/3xj)ds/ V 

where  V^UiUjlg  is  the  initial  mean  kinetic  energy  at  the  point  of  release  of  the  fluid  parti¬ 
cle,  ds  is  the  infinitesimal  arclength  along  the  streamline  and  V  =  V  ujuj.  The  separate 
contributions  of  the  integrands  -(ui/p)3p/3xj  and  (uj/p)3Tjj/3xj  provide  a  quantitative 
measure  of  the  inviscid  and  viscous  contributions  to  the  rate  of  change  of  the  mean 
kinetic  energy. 


Determination  of  Flowfield  Models 


During  the  past  year,  a  significant  effort  was  focused  on  the  determination  of  the  flowfield 
structure  for  the  3-D  swept  compression  using  the  computed  flowfields  obtained  with  the 
Baldwin-Lomax  turbulence  model.  These  flowfields  had  been  validated  by  comparison  with  the 
experimental  data  of  Ketchum  and  Bogdonoff  (Knight  et  al  1988),  Settles  and  McKenzie  (Set¬ 
tles  et  al  1986),  Settles  and  Teng  (Settles  and  Teng  1984),  and  Ruderich  and  Mao  (Knight  et  al 
1987b).  Details  of  the  validation  are  presented  in  Knight  et  al  (1987b)  and  Knight  et  al  (1988). 
The  flowfield  analysis  was  performed  using  the  software  and  graphical  techniques  described 
above.  Efforts  were  focused  on  the  3-D  swept  compression  comer  for  tips  (a,  X.)  =  (24  ,  40) 
and  (24,  60)  deg  configurations  at  Mach  3  for  Reg^  =  1.4xl05  to  9x10  .  The  principal  con¬ 
clusions,  described  in  detail  in  Knight  et  al  (1988),  are  the  following  : 


1.  Vortical  Structure 


jCLa, 


The  3-D  swept  compression  corner  flowfield  is  dominated  by  a  large  vortical  structure  for 
both  the  (a,X)  =  (24,  40)  and  (24,  60)  deg  configurations  at  Mach  3.  The  overall  struc¬ 
ture  is  characterized  by  two  surfaces  (Fig.  3),  namely,  1)  the  surface  emanating  from  the 
line  of  separation  or  coalescence  (denoted  as  ’Surface  No.  1’),  and  2)  the  surface  inter¬ 
secting  the  line  of  attachment  or  divergence  (denoted  as  ’Surface  No.  2’).  The  fluid  con¬ 
tained  between  the  wall  and  the  second  surface  is  entrained  within  the  vortex,  while  the 
fluid  above  the  second  surface  passes  over  the  vortex  and  moves  up  the  ramp. 


2.  Viscous  and  Inviscid  Effects 


The  contributions  to  the  rate  of  change  of  the  mean  kinetic  energy  were  evaluated  for  a 


series  of  streamlines  for  the  (a,  X)  =  (24,  40)  deg  configuration  at  Reg^  =  8. lx 


10“ 
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detailed  description  of  this  analysis  is  presented  in  Knight  et  al  (1988).  It  was  observed 
that  the  rate  of  change  of  the  mean  kinetic  energy  was  principally  due  to  inviscid  effects 
except  for  streamlines  originating  very  dose  to  the  surface  in  the  upstream  boundary  layer 
(i.e.,  y  <  0.058^).  A  similar  dialysis  is  in  progress  for  the  (a,  X)  =  (24,  40)  deg 
configuration  for  Ret  =  2.6xl05  and  the  (a,  X)  =  (24  ,  60)  deg  configuration  at  Reg^ 
=  1.4x10  and  9x10  . 

These  results  are  consistent  with  the  previous  observation  regarding  the  insensitivity 
of  the  computed  pitot  pressure  and  yaw  angle  profiles  to  the  turbulence  model  employed, 
except  in  the  immediate  vicinity  of  the  surface.  The  3-D  turbulent  interactions  investi¬ 
gated  are  therefore  principally  rotational  and  inviscid  within  the  boundary  layer,  except 
close  to  the  wall. 

Research  Task  No.  3:  Investigate  Methods  for  Control  and 

Modification  of  3-D  Turbulent  Interactions 

The  principal  function  of  a  conventional  high  speed  aircraft  inlet  is  to  provide  uniform, 
high  total  pressure,  subsonic  flow  at  the  compressor  face.  It  has  long  been  observed  that  boun¬ 
dary  layer  separation  within  the  inlet  can  lead  to  regions  of  low  total  pressure  recovery  at  the 
inlet  surface,  and  subsequent  degradation  of  engine  performance.  Boundary  layer  bleed  has 
traditionally  been  employed  in  high  speed  aircraft  inlets  to  prevent  boundary  layer  separation. 

The  development  of  accurate  theoretical  models  for  the  3-D  sharp  fin  and  swept 
compression  corner  configurations  provides  the  opportunity  for  investigation  of  methods  to 
control  and  modify  these  interactions.  During  FY87  efforts  have  focused  on  the  3-D  sharp  fin 
configuration  and  modification  by  surface  bleed. 

The  focus  of  the  activity  on  modification  by  surface  bleed  is  the  3-D  sharp  fin  interac¬ 
tion.  This  configuration  is  analogous  to  the  sidewall  interaction  in  a  high  speed  inle»,  with  the 
sharp  fin  representing  the  ramp  surface,  and  the  flat  plate  representing  the  inlet  sidewall.  It  is 
known  that  the  effect  of  the  sidewall  interaction  is  to  generate  substantial  spanwise  flow  which, 
upon  interaction  with  the  cowl  (opposite  the  ramp),  causes  substantial  low  speed  flow  within 
the  cowl  boundary  layer  and  subsequent  low  pressure  recovery  at  the  compressor  face.  The 
objective  of  this  research  activity,  therefore,  is  to  determine  the  effect  of  surface  bleed  on  the 
flat  plate  (inlet  sidewall)  on  the  large  vortical  structure. 

A  series  of  computations  have  been  performed  for  the  3-D  sharp  fin  at  Mach  3  for  a  fin 
angle  of  20  deg  at  Reg^  =  9x10  .  Boundary  layer  bleed  was  provided  in  a  triangular  region 
between  the  line  of  upstream  influence  (as  defined  by  the  surface  pressure)  and  the  theoretical 
inviscid  shock  location.  Three  separate  bleed  flow  rates  (1%,  2.5%  and  5%  based  on  p^U^) 
were  employed.  These  represent  moderate  to  high  bleed  flow  rates  on  basis  of  current  design 
practice. 


A  detailed  description  of  the  results  is  provided  in  Gaitonde  and  Knight  (1988)*.  The 
principal  conclusions  are  : 


Included  m  Section  VI 


1.  Effect  on  Vortical  Structure 


The  bleed  has  a  negligible  effect  on  the  vortical  structure.  The  large  vortex,  representing 
the  principal  structure  of  the  3-D  turbulent  interaction,  is  essentially  unaffected  in  size 
and  strength.  This  represents  a  significant,  and  indeed  surprising,  result.  In  particular,  it 
brings  into  question  the  efficacy  of  boundary  layer  bleed  for  sidewall  interactions  in  high 
speed  inlets. 

2.  Effect  on  Compression  Wave  Ahead  of  Shock 

The  compression  wave  system,  located  upstream  of  the  shock,  is  tightened  by  the  surface 
bleed,  which  effectively  reduces  the  displacement  thickness  of  the  interaction.  The  reduc¬ 
tion  in  size  of  the  compression  region  is  evident  in  the  pitot  pressure  and  static  pressure 
profiles. 

3.  Effect  on  Surface  Skin  Friction  Lines 

The  line  of  coalescence  is  observed  to  move  towards  the  downstream  boundary  of  the 
bleed  region.  The  surface  skin  friction  lines  between  the  lines  of  coalescence  and  diver¬ 
gence  show  the  same  qualitative  behavior  as  observed  for  the  no  bleed  case. 

Research  Task  No.  4:  Develop  a  Unified  Understanding  of  the 

Flowfield  Structure  of  3-D  Turbulent  Interactions 

On  the  basis  of  detailed  examination  of  particle  traces  for  the  3-D  swept  compression 
corner  at  (a,  X)  =  (24,  40)  and  (24,  60)  deg  during  the  present  year,  it  is  evident  that  the 
flowfield  structure  for  these  3-D  swept  compression  comer  configurations  and  the  3-D  sharp  fin 
are  qualitatively  identical.  This  represents  the  first  step  in  development  of  a  more  unified  theory 
of  flowfield  structures  for  3-D  turbulent  interactions. 


B.  Research  Program  for  FY88  (1  October  1987  -  30  September  1988) 

The  research  program  for  FY88  is  characterized  by  two  major  themes,  namely,  i)  the 
continuation  of  the  research  effort  in  supersonic  shock  wave-turbulent  boundary  layer  interac¬ 
tion,  and  ii)  the  initiation  of  a  significant  research  activity  in  3-D  hypersonic  shock  wave- 
turbulent  boundary  layer  interaction.  The  basic  four  objectives,  outlined  in  the  original  propo¬ 
sal  (Knight  1986b)  and  described  in  Section  I,  are  pertinent  to  both  the  supersonic  and  hyper¬ 
sonic  research  activity.  The  research  tasks  for  the  remainder  of  FY88  are  described  below  in 
terms  of  the  four  objectives.  It  is  noted  that  some  of  the  activities  are  expected  to  continue  into 
FY89. 


Research  Task  No.  1  :  Develop  and  Validate  Theoretical  Models 

for  3-D  Shock  Wave-Turbulent  Boundary  Layer  Interactions 

The  recent  establishment  of  the  National  Aerospace  Plane  project  (NASP)  has  led  to 
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renewed  interest  in  hypersonic  flight.  The  phenomena  of  three  dimensional  shock  wave- 
turbulent  boundary  layer  interaction  is  prevalent  in  hypersonic  flight,  occurring  in  a  variety  of 
practical  vehicle  configurations.  Despite  an  active  research  effort  in  hypersonic  flows  in  the  late 
sixties  and  early  seventies,  a  complete  understanding  of  3-D  hypersonic  turbulent  interactions 
has  not  been  achieved  (Holden  1986).  Recently,  Horstman  (1987a)  has  studied  2-D  hyper¬ 
sonic  turbulent  boundary  layer  interactions  using  the  k-e  turbulence  model  and  three  different 
modifications  to  account  for  hypersonic  compressibility  effects.  He  concluded  that  no  single 
turbulence  model  is  capable  of  accurate  prediction  of  the  experimental  data.  The  applicability 
of  this  result  to  3-D  hypersonic  turbulent  interactions,  however,  is  not  clear.  Recent  experi¬ 
mental  and  theoretical  results  at  Mach  3  have  shown  that  simple  turbulence  models  (i.e.,  zero 
equation  and  two  equation  eddy  viscosity  models)  are  capable  of  accurate  simulation  of  a 
variety  of  3-D  turbulent  interactions  (Knight  et  al  1987a,  Knight  et  al  1987b,  Knight  et  al 
1988),  while  failing  to  accurately  reproduce  the  limiting  2-D  turbulent  interaction  (Ong  and 
Knight,  1986).  The  success  in  3-D  supersonic  turbulent  interactions  was  attributed  to  the 
nature  of  the  interaction,  specifically,  the  flowfields  investigated  were  observed  to  be  essentially 
inviscid  and  rotational,  except  in  a  small  fraction  of  the  turbulent  boundary  layer  adjacent  to  the 
surface  where  the  effects  of  turbulence  were  observed  to  be  significant.  A  collaborative  experi¬ 
mental  and  theoretical  investigation  of  3-D  hypersonic  turbulent  interactions  is  therefore 
required. 


The  proposed  research  project  involves  the  numerical  simulation  of  the  3-D  hypersonic 
turbulent  interaction  generated  by  a  sharp  fin.  The  flow  configuration  is  shown  in  Fig.  1.  An 
equilibrium  turbulent  boundary  layer  develops  on  the  flat  plate.  The  deflection  of  the  fin 
creates  an  oblique  shock  wave,  which  interacts  with  the  turbulent  boundary  layer  on  the  flat 
plate.  The  selection  of  this  configuration  was  motivated  by  several  factors.  First,  the  recent 
collaborative  experimental  and  theoretical  effort  involving  the  principal  investigator,  NASA 
Ames  Research  Center  and  Princeton  Gas  Dynamics  Laboratory  (Knight  et  al  1987a)  has 
developed  a  detailed  understanding  of  the  fiowfield  structure  of  the  3-D  sharp  fin  turbulent 
interaction  at  Mach  3.  Second,  a  series  of  3-D  sharp  fin  experiments  are  planned  at  hypersonic 
speeds  (Mach  7  and  10)  at  NASA  Ames  Research  Center  in  FY88  and  FY89  (Horstman 
1987b). 


A  series  of  computations  will  be  performed  for  the  3-D  sharp  fin  configuration  at  Mach 
7  at  fixed  Reynolds  number  for  several  fin  angles  a.  Additional  computations  may  be  per¬ 
formed  at  Mach  10.  The  flow  conditions  will  be  selected  to  agree  with  the  planned  experiments 
at  NASA  Ames.  Although  the  precise  requirements  for  the  number  of  grid  points  and  cpu 
time  depends  on  the  specific  case,  experience  has  indicated  that  typically  100,000  grid  points 
and  10  CYBER  205  cpu  hrs  are  required  (Knight  et  al  1987b)  for  each  case. 


The  computed  results  will  be  compared  with  the  experimental  data.  Present  plans 
include  measurements  of  surface  pressure  and  heat  transfer,  surface  oil  film  visualization,  and 
boundary  layer  profiles  of  pitot  and  static  pressure,  velocity,  and  turbulence  fluctuations.  The 
data  and  graphics  postprocessing  will  be  performed  at  the  Rutgers  College  of  Engineering 
Supercomputer  Remote  Access  and  Graphics  Facility. 
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Research  Task  No.  2: 


Determination  of  Physical  Structure 

of  3-D  Turbulent  Interactions  for  Specific  Geometries 


Recent  collaborative  efforts  involving  the  principal  investigator,  NASA  Ames  Research 
Center  and  the  Princeton  Gas  Dynamics  Laboratory  have  focused  on  the  understanding  of  the 
flowfield  structure  for  two  specific  dimensionless  configurations,  namely,  the  sharp  fin  (Fig.  1) 
and  the  swept  compression  comer  (Fig.  2).  On  the  basis  of  these  investigations,  conducted  at 
Mach  3  over  a  Reynolds  number  range  from  2.3xl(T  to  9x10  ,  a  view  of  the  overall  flowfield 
structure  has  emerged.  The  flowfield  is  inviscid  and  rotational  over  a  majority  of  the  interaction 
domain,  except  within  a  narrow  region  adjacent  to  the  surface  where  the  effects  of  turbulent 
mixing  are  significant.  For  both  the  sharp  fin  at  Mach  3  (and  a  up  to  20  deg)  and  the  swept 
compression  comer  at  Mach  3  (and  (a.A.)  =  (24,40)  and  (24,60)  deg),  the  flowfield  is  dom¬ 
inated  by  a  large  vortical  structure. 

A  series  of  experiments  were  performed  recently  at  the  Princeton  Gas  Dynamics 
Laboratory  to  investigate  the  effect  of  the  shock  generating  geometry  on  the  interaction  (Kim- 
mel  1987).  Three  configurations  were  selected,  namely,  the  sharp  fin  (Fig.  1),  swept  compres¬ 
sion  comer  (Fig.  2),  and  semicone  (Fig.  4).  The  geometrical  characteristics  were  selected  to 
obtain  similar  shock  strengths.  Specifically,  the  sharp  fin  angle  a  =  17.5  deg,  the  angle  of 
compression  and  sweep  for  the  swept  compression  comer  areas  30  deg  and  1*  60  deg,  and 
the  semicone  half  angle  v  ■  25  det respectively.  Each  configuration  was  examined  at  Rey¬ 
nolds  numbers  of  1.6xlOyand  9-lxlO5.  For  the  tower  Reynolds  number,  detailed  surface  pres¬ 
sure  and  kerosene-lampblack  flow  visualization  was  performed.  At  the  higher  Reynolds 
number,  experimental  data  include  both  surface  pressure  and  kerosene-lampblack  flow  visuali¬ 
zation,  and  boundary  layer  profiles  of  pitot  pressure,  static  pressure,  and  yaw  angle. 

The  three  configurations  exhibited  conical  similarity,  with  the  flow  surface  and  flowfield 
variables  exhibiting  two  distinct  regions.  The  first  region,  extending  from  upstream  to  a  point 
downstream  of  the  shock  system,  displayed  a  ’quasi-con ical  free  interaction'.  In  this  region,  for 
example,  the  experimental  pressure  distribution,  displayed  along  a  line  normal  to  the  inviscid 
shock,  was  virtually  identical  for  the  three  configurations  despite  their  significant  geometrical 
differences.  A  second  region,  downstream  of  the  first  region  and  denoted  the  'model  dominated 
region’,  showed  no  similarity  of  profiles. 

The  explanation  for  the  quasi-conical  free  interaction  is  not  clear.  It  is  anticipated  that 
all  three  flowfields  display  a  large  vortical  structure.  The  sharp  fin  at  a  *  17.5  deg  is  inter¬ 
mediate  to  a  series  of  computed  flows  (a  =  10  and  20  deg)  which  have  displayed  a  large  vorti¬ 
cal  structure  (Knight  et  al  1987a).  The  swept  compression  comer  at  (a,X)  *  (30,60)  deg  is 
close  to  the  previous  computation  of  Knight  et  al  at  (a,\)  =  (24,60)  deg.  The  semicone,  how¬ 
ever,  has  not  been  computed  previously.  It  is  unclear,  however,  why  this  vortical  structure 
should  yield  nearly  identical  pressure  distributions  in  the  quasi-conical  free  interaction  region. 

The  proposed  research  will  focus  on  the  three  configurations,  namely,  the  sharp  fin  at  a 
=  17.5  deg,  the  semicone  at  y  =  25  deg,  and  the  swept  compression  comer  at  (a,\)  = 
(30,601  deg.  Computations  will  be  performed  at  Mach  3  for  Reynolds  numbers  of  1.6xl03  and 
9.1x10  .  Although  the  precise  requirements  for  the  number  of  grid  points  and  cpu  time 
depends  on  the  specific  case,  experience  has  indicated  that  typically  100,000  grid  points  and  10 
CYBER  205  cpu  hrs  are  required  (Knight  et  al  1987b)  for  each  case. 

The  computed  results  will  be  compared  with  the  experimental  data  of  Kimmei  (1987). 


-13- 


The  computed  flowfields  will  be  analyzed  to  determine  the  flowfield  structure,  and  examine  the 
reasons  for  the  quasi-conical  similarity.  The  data  and  graphics  postprocessing  will  be  performed 
at  the  Rutgers  College  of  Engineering  Supercomputer  Remote  Access  and  Graphics  Facility. 
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Research  Task  No.  3: 


Investigate  Methods  for  Control  and 
Modification  of  3-D  Turbulent  Interactions 


The  research  activity  on  control  and  modification  of  3-D  turbulent  interactions  is  con¬ 
tinued  in  FY88.  Two  specific  activities  are  in  progress  : 

1.  3-D  Sharp  Fin  in  the  Presence  of  Bleed 

The  research  effort  on  the  effects  of  boundary  layer  bleed  for  the  3-D  sharp  fin 
configuration  focused  on  a  specific  bleed  configuration  in  FY87,  namely,  a  triangular  bleed 
region  confined  between  the  line  of  upstream  influence  and  the  inviscid  shock  (Section 
II. A).  A  series  of  computations  were  performed  at  Mach  3  for  a  fin  angle  of  20  deg  at 
Re*  =  9x10s.  During  the  present  year,  the  same  configuration  will  be  examined  with  a 


different  bleed  schedule,  namely,  a  bleed  region  confined  between  the  inviscid  shock  and 
the  fin  surface.  These  two  configurations  will  therefore  provide  an  examination  of  the 
effects  of  bleed  throughout  the  entire  interaction  region. 

3-D  Intersecting  Shock  Wave-Turbulent  Boundary  Layer  Interaction 

A  novel  concept  was  proposed  by  Mee  and  Stalker  (1987)  to  control  3-D  shock  wave- 
turbulent  boundary  layer  interactions.  They  argued  that  "intersecting  shock  interactions 
can  produce  a  given  overall  pressure  rise  with  less  likelihood  of  separation  than  an 
equivalent  strength  single  shock  interaction".  A  model  of  the  flow  is  displayed  in  Fig.  5. 
Two  sharp  fins,  mounted  normal  to  the  flat  plate,  generate  oblique  shock  waves  which 
intersect.  The  fin  angles  may  be  equal  or  dissimilar,  yielding  a  symmetric  or  asymmetric 
interaction,  respective.  The  initial  3-D  shock  wave-turbulent  boundary  layer  interaction 
associated  with  each  oblique  shock  generates  a  vortical  structure  whose  strength  is  depen¬ 
dent  on  the  flow  conditions  and  fin  angle  (i.e.,  initial  pressure  rise).  The  sense  of  rotation 
of  the  flow  within  the  vortical  structures  are  of  opposite  sign. 

The  concept  of  Mee  and  Stalker  is  unproven,  however,  and  additional  computational  and 
experimental  research  is  required  to  determine  its  validity.  In  FY  88,  two  specific  compu¬ 
tations  are  planned.  First,  a  calculation  will  be  performed*  for  the  strongest  symmetric 
intersecting  shock  configuration  of  Mee  and  Stalker.  This  case  corresponds  to  an  incom¬ 
ing  Mach  number  of  1.85,  Reynolds  number  Reg^  =  7.8x10,  and  symmetric  sharp  fins 
with  an  angle  of  5  deg.  This  represents  a  rather  weak  interaction,  with  an  overall  pressure 
rise  across  the  intersection  shocks  Pgnai/PUDStream  e9ual  t0  1-7-  Second,  a  series  of  com¬ 
putations  will  be  performed  for  the  intersecting  shock  configurations  currently  under  study 
at  the  Princeton  Gas  Dynamics  Laboratory.  These  <jases  correspond  to  an  incoming  Mach 
number  of  1.85,  Reynolds  number  Reg^  =  1x10s  to  2x10s,  and  symmetric  sharp  fins 
with  angles  of  4,  6  and  8  deg.  The  experimental  program  at  the  Princeton  Gas  Dynamics 
Laboratory  includes  measurements  of  surface  pressure  and  surface  flow  visualization. 
These  configurations  represent  significantly  stronger  interactions,  with  overall  pressure 


Pfinal^P 


upstream 


=  1.8,  2.4,  3.1  for  fin  angles  of  4,  6,  and  8  deg. 


*  This  computation  has  been  completed,  and  is  included  in  the  paper  by  Gaitonde  and  Knight  in  Section 


Research  Task  No.  4:  Develop  a  Unified  Understanding  of  the 

Flowfield  Structure  of  3-D  Turbulent  Interactions 

The  research  program  continues  to  focus  on  the  objective  of  developing  a  unified 
understanding  of  the  flowfield  structure  of  3-D  turbulent  interactions.  In  the  previous  year,  it 
was  observed  that  the  flowfield  structure  for  specific  configurations  of  the  3-D  sharp  fin  and  3-D 
swept  compression  corner  interactions  was  similar  (i.e.,  a  large  vortical  structure).  Efforts  con¬ 
tinue  to  identify  similarities  in  the  different  configurations  examined  under  the  previous  three 
objectives. 
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Section  III.  Publications  and  Scientific  Interactions 


Written  Publications  -  Cumulative  Chronological  List 

1  October  1981  -  30  September  1982 

a.  Knight,  D.,  "Application  of  Curvilinear  Coordinate  Generation  Techniques  to  the 
Computation  of  Internal  Flows",  in  Numerical  Grid  Generation  -  Proceedings  of  a 
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Fig.  5  3-D  Intersecting  Shock  Waves 
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Abstract 

A  combined  experimental  and  theoretical 
investigation  has  continued  for  the  3-D  shock 
wave-turbulent  boundary  layer  interaction  gen¬ 
erated  by  a  swept  compression  comer.  The 
compression  comer  geometry  is  characterized  by 
the  angle  of  streamwise  compression  a  and  the 
angle  of  sweep  X.  The  present  study  focuses  on 
the  (a.X)  =  (24,40)  deg  configuration  at  Reg^ 

*  2.6x10s  and  8.1x10s.  Recent  computations 
have  been  performed  for  this  configuration  using 
the  Baldwin -Lomax  algebraic  turbulent  eddy 
viscosity  model.  Previous  computations  were 
performed  using  the  Cebeci-Smith  algebraic  tur¬ 
bulent  eddy  viscosity  model  and  the  Jones- 
Launder  two-equation  (k-e)  model.  The  results 
using  these  three  models  are  compared  in  the 
present  paper  with  the  experimental  measure¬ 
ments  of  Settles  and  McKenzie  for  (a,X)  « 
(24,40)  deg  at  Reg_  -  2.6x10s  and  8.1x10s, 
and  the  experimental  measurements  of  Ketchum 
and  Bogdonoff  at  Reg..  «  8.1x10s.  The  com¬ 
puted  flowfields  are  generally  in  good  agreement 
with  experimental  profiles  of  pitot  pressure  and 
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yaw  angle  at  both  Reynolds  numbers,  and  are 
relatively  insensitive  to  the  turbulence  model 
employed.  The  computed  surface  pressure 
underpredicts  the  upstream  influence  of  the 
interaction.  The  flowfield,  examined  using  parti¬ 
cle  tracing,  manifests  a  large  vortical  structure. 
The  quantitative  details  of  the  flowfield  structure 
are  examined  for  the  (a.X)  *  (24,40)  deg 

configuration,  and  compared  to  previous  results 
for  the  (a.X)  *  (24,60)  deg  configuration.  The 
rate  of  change  of  the  mean  kinetic  energy  along 
a  streamline  is  investigated.  The  principal  contri¬ 
bution  is  attributable  to  inviscid  effects  except 
for  streamlines  originating  very  close  to  the  wall. 
This  examination  provides  direct  evidence  that 
the  flowfield  structure  is  rotational  and  inviscid 
except  in  the  immediate  vicinity  of  the  surface. 

Introduction 

The  interaction  between  shock  waves  and 
turbulent  boundary  layers  is  an  important  prob¬ 
lem  in  modem  fluid  mechanics,  with  applications 
to  high  speed  aerodynamics  and  propulsion.  The 
earliest  observations  were  apparently  made  by 
Ferri^,  and  the  first  systematic  investigations 
were  performed  by  Liepmanns  and  Ackeret4.  In 
a  review  of  3-D  supersonic  turbulent  interactions 
in  1969,  Greene  remarkeds, 

"We  have  seen  that  present  understanding  of 
shock  and  boundary  layer  interactions  in  two- 
dimensional  flow  is  not  entirely  satisfactory.  In 
general  three-dimensional  flows  it  scarcely  exists 
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in  any  quantitative  sense." 

Within  the  past  fifteen  years,  a  significant  effort 
has  been  focused  on  3-D  shock  wave-turbulent 
boundary  layer  interactions  (denoted  "3-D  tur¬ 
bulent  interactions").  The  research  has 
employed  experimental,  computational,  and 
analytical  techniques.  A  general  review  of  swept 
shock  wave-turbulent  boundary  layer  interactions 
has  been  presented  by  Settles  and  Dolling^. 

A  series  of  collaborative  experimental  and 
theoretical  investigations  by  the  authors  has 
focused  on  two  specific  dimensionless  geometries 
for  3-D  shock  wave-turbulent  boundary  layer 
interactions,  namely,  the  3-D  sharp  fin  (Fig.  1) 
and  the  3-D  swept  compression  comer  (Fig.  2). 
The  theoretical  model  is  the  Reynolds-averaged 
compressible  Navier-Stokes  equations,  with  tur¬ 
bulence  incorporated  using  an  turbulent  eddy 
viscosity.  Several  different  turbulent  eddy 
viscosity  models  are  employed,  including  the 
Bald  win -Lomax  ^  and  Cebeci-Smith"^  algebraic 
models  and  the  Jones-Launder®  or  k-e  two  equa¬ 
tion  model.  The  governing  equations  are  solved 
using  two  different  numerical  algorithms.  The 
experiments  utilize  the  high  Reynolds  number 
wind  tunnel  at  the  Princeton  Gas  Dynamics 
Laboratory.  The  3-D  sharp  fin  has  been 
studied** at  Mach  3  for  Reynolds  numbers 
Re6—  ”  2.5x10s  to  9x10s  and  fin  angle  a  up  to 
20  deg.  It  was  observed  that  the  computed 
results  were  in  good  agreement  with  the  experi¬ 
mental  data  for  surface  pressure,  and  boundary 
layer  profiles  of  pitot  pressure,  static  pressure 
and  yaw  angle.  The  computed  flowfields  were 
found  to  be  insensitive  to  the  turbulence  model 
employed,  except  in  a  small  fraction  of  the  boun¬ 
dary  layer  adjacent  to  the  surface,  where  the 
computed  profiles  displayed  modest  differences 
amongst  the  predictions  of  the  various  tur¬ 
bulence  models  and  the  experiment.  It  was  con¬ 
cluded,  therefore,  that  the  3-D  sharp  fin  interac¬ 
tion  was  principally  in  viscid  and  rotational, 
except  in  the  aforementioned  narrow  region  adja¬ 
cent  to  the  surface.  The  flowfield  was  observed 
to  be  dominated  by  a  large  vortical  structure. 

In  a  previous  paper*4,  the  3-D  swept 
compression  comer  for  the  (a,X)  «  (24,60) 

configuration  was  examined  at  Rej^  =  1.4x10s 
to  9x10s.  Computations  were  performed  using 
the  Baldwin -Lorn  ax  and  Jones-Launder  tur¬ 
bulence  models,  and  the  results  compared  with 


the  experimental  data  of  Ruderich,  Mao  and 
Bogdonoff.  The  computed  flowfields  for  Re^ 
-  9X105  were  found  to  be  in  good  agreement 
with  the  experimental  data  for  surface  pressure, 
and  boundary  layer  profiles  of  pitot  pressure  and 
yaw  angle.  The  computed  surface  pressure  for 
the  lower  Reynolds  number  case,  however, 
displayed  significant  (i.e.,  25%)  differences  with 
experiment.  The  computed  flowfields  for  the 
-  9x10s  case  displayed  a  large  vortical 
structure,  with  a  topology  similar  to  the  3-D 
sharp  fin. 

The  present  paper  focuses  on  the  3-D 
swept  compression  comer  for  the  (a,X)  = 
(24,40)  at  Re^  *  2.6x10s  and  8.1x10s. 

Recent  computations  have  been  performed  for 
this  configuration  using  the  Baldwin-Lomax 
model.  Previous  calculations  were  performed  by 
Horstman*S  using  the  Cebeci-Smith  and  Jones- 
Launder  turbulence  models.  The  objectives  of 
the  present  effort  are  : 

1.  Accuracy  of  the  theoretical  models 

The  efficacy  of  the  Baldwin-Lomax  model  for  the 
(a,X)  =  (24,40)  deg  configuration  is  examined 
through  comparison  with  the  present  experimen¬ 
tal  data  of  Ketchum  and  Bogdonoff  and  the  prior 
experimental  study  of  Settles  et  aJls  for  surface 
pressure,  and  boundary  layer  profiles  of  pitot 
pressure  and  yaw  angle.  These  results  are  also 
compared  with  the  previous  computations  of 
Horstman  using  the  Cebeci-Smith  and  Jones- 
Launder  models. 

2.  Determination  of  the  Flowfield  Structure 

Provided  reasonable  agreement  is  obtained 
between  the  computed  and  measured  data,  the 
computed  flowfields  can  be  employed  to  ascertain 
the  basic  flowfield  structure.  The  (a,X)  = 
(24,40)  deg  configuration  is  examined,  with 
emphasis  on  quantifying  the  nature  of  the  topol¬ 
ogy  of  the  flow.  This  includes  the  description  of 
the  distinct  regions  of  the  flowfield  and  an  exam¬ 
ination  of  the  viscous  and  inviscid  contributions 
to  the  rate  of  change  of  the  mean  kinetic  energy. 

Description  of  Experiments 

The  experiments  were  conducted  in  the 
high  Reynolds  number  wind  tunnel  at  the 
Princeton  Gas  Dynamics  Laboratory.  A  descrip¬ 
tion  of  the  facility  is  provided  in  Vas  and  Bog 
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donoff1^.  The  facility  is  a  blowdown  tunnel  with 
a  20  cm  x  20  cm  test  section  and  a  nominal 
Mach  3  test  section.  The  total  pressure  for  the 
(a.X)  =  (24,40)  deg  configuration  was  683  kPa 
t.  1.84b  and  690  kPa  -tl-39b  for  Reg^.  = 
2.6x10^  and  8.1x10^,  respectively.  The  total 
temperature  for  the  (a.X)  =  (24,40)  deg 
configuration  was  265.4  deg  K  -tl-^fc  for  both 
Reynolds  numbers.  The  surface  temperature 
was  near  adiabatic.  The  incoming  boundary  layer 
was  an  equilibrium,  two-dimensional  turbulent 
boundary  layer  which  has  been  extensively  sur¬ 
veyed^-  *7  and  observed  to  satisfy  the  Law  of 
the  Wall  and  Wake  18,19.  The  nominal  boundary 
layer  thickness  6_  =  0.42  cm  and  1.27  cm  for 
Reg^  =  2.6x10^  and  8.1x10^,  respectively. 

Detailed  boundary  layer  profile  measurements 
were  performed  for  the  5—  =  127  cm 
configuration.  The  size  of  the  instrumentation 
limited  the  extent  of  the  boundary  layer  profile 
data  for  the  5  *  0.42  cm  case. 


Description  of  Computations 


A  total  of  three  different  theoretical 
models  are  employed  for  the  3-D  swept  compres¬ 
sion  comer  at  (a.X)  «  (24,40)  deg.  The 

governing  equations  are  the  full  mean  three- 
dimensional  compressible  Navier-Stokes  equa¬ 
tions  using  mass-averaged  variables70  and  strong 
conservation  form71.  Three  separate  turbulent 
eddy  viscosity  models  are  employed,  namely,  the 
algebraic  models  of  Baldwin-Lomax**  and 
Cebeci-Smith7,  and  the  two-equation  (k-E) 
model  of  Jones-Launder**.  The  molecular 
dynamic  viscosity  is  specified  by  Sutherland's 
law.  The  molecular  and  turbulent  Prandtl 
numbers  are  0.73  and  0.9,  respectively. 


The  computational  domain  is  shown  in 
Fig.  3.  An  equilibrium  turbulent  boundary  layer 
profile  is  prescribed  at  the  upstream  boundary 
ABHG  satisfying  the  measured  momentum 
thickness.  It  is  noted  that  the  upstream  boun¬ 
dary  is  parallel  to  the  comer  line,  and  hence 
swept  downstream  at  an  angle  X  relative  to  the 
spanwise  coordinate  z.  The  calculations  employ¬ 
ing  the  Baldwin-Lomax  model  incorporated  the 
variation  in  upstream  boundary  layer  thickness 
on  ABHG  due  to  the  sweep  of  ABHG.  The  cal¬ 
culations  using  the  Cebeci-Smith  and  Jones- 
Launder  models  employed  an  average  boundary 
layer  profile.  Previous  computations  at  (a.X)  * 
(24,60)  deg  have  shown  that  the  3-D  interaction 
is  insensitive  to  moderate  changes  in  the 


upstream  boundary  layer  thickness1*.  At  the 
plane  of  symmetry  AFEKLG,  the  normal  com¬ 
ponent  of  the  velocity  is  set  to  zero,  and  the  nor¬ 
mal  gradient  of  the  remaining  flow  variables  is 
zero.  On  the  solid  surface  ABCDEF,  the  velo¬ 
city  is  zero,  a  fixed  surface  temperature  (near 
adiabatic  conditions)  is  prescribed,  and  the  nor¬ 
mal  derivative  of  the  static  pressure  is  set  to 
zero.  At  the  right  boundary  BCDJIH,  the  gra¬ 
dient  of  the  flow  variables  in  a  direction  parallel 
to  the  comer  line  is  set  to  zero.  On  the  down¬ 
stream  boundary  EDJK,  the  gradient  of  the  flow 
variables  is  set  to  zero  along  a  direction  parallel 
to  the  x-axis  and  compression  ramp  surface. 
The  upper  boundary  GHIJKL  is  fixed  sufficiently 
far  from  the  interaction  that  the  flowfield  is 
undisturbed  and  uniform  freestream  conditions 
are  applied. 


The  computations  using  the  Cebeci-Smith 
and  Jones-Launder  models  were  solved  using  the 
hybrid  explicit-implicit  algorithm  of  MacCor- 
mack77.  The  computations  using  the  Baldwin- 
Lomax  model  were  solved  using  the  hybrid 
explicit-implicit  algorithm  of  Knight10.  Both 
algorithms  have  been  highly  vectorized  for  exe¬ 
cution  on  CRAY  and  CYBER  205  computers, 
respectively. 


The  results  of  the  three  different  theoreti¬ 
cal  models  for  the  (a.X)  =  (24.40)  deg  compres¬ 
sion  comer  at  Rej^  =  2.6x10^  are  compared 
with  the  experimental  data  of  Settles  and 
McKenzie.  The  computations  employing  the 
Cebeci-Smith  and  Jones-Launder  models  utilized 
a  grid  of  27,200  points.  The  x-grid  spacing  was 
concentrated  near  the  comer,  with  a  minimum 
Ax/5—  =  0.18,  and  a  maximum  Ax/5.„  *  1.2 
at  the  upstream  and  downstream  boundaries. 
The  y-grid  spacing  was  concentrated  near  the 
surface  with  the  minimum  Ay/8„  =  2.78x10'*, 
and  the  maximum  Ay/6..  *  1.39  at  the  upper 
surface  GHIJKL.  The  grid  spacing  adjacent  to 
the  surface  satisfied  Ayu«/vw  <  1.0.  A  uniform 
z-grid  spacing  Az/5.,  =  0.727  was  employed. 
The  computation  employing  the  Baldwin-Lomax 
model  utilized  a  grid  of  244,400  points.  The  x- 
grid  spacing  was  concentrated  near  the  corner, 
with  a  minimum  Ax/5..  *  0.157.  The  max¬ 
imum  Ax/8..  =  0-5 1 5  and  0.788  at  the 

upstream  and  downstream  boundaries,  respec¬ 
tively.  The  y-grid  spacing  was  concentrated  near 
the  surface  with  the  minimum  Ay/ 5..  = 

1.03x10'^,  and  the  maximum  Ay/8.,  =  0.684  at 


the  upper  surface  GHIJKL.  The  z-grid  spacing 
was  concentrated  near  the  symmetry  plane 
AFEKLG,  with  a  minimum  Az/fia.  =  1.82x10" 

2  and  a  maximum  Az/6_  =  0.694  near  the 
boundary  BCDJIH. 

The  results  of  the  three  different  theoreti¬ 
cal  models  for  the  (a,X)  =  (24,40)  deg  compres¬ 
sion  corner  at  Res*,  *  8.1x10^  are  compared 
with  the  experimental  data  of  i)  Settles  and 
McKenzie,  and  ii)  Ketchum  and  Bogdonoff.  The 
computations  employing  the  Cebeci-Smith  and 
Jones-Launder  models  utilized  a  grid  of  28,000 
points.  The  x-grid  spacing  was  concentrated  near 
the  comer,  with  a  minimum  Ax/5*,  =  0.3,  and 
a  maximum  Ax/5^  =  1.0  at  the  upstream  and 
downstream  boundaries.  The  y-grid  spacing  was 
concentrated  near  the  surface  with  the  minimum 
Ay/8*,  *  1.02xl0"4,  and  the  maximum  Ay/8.. 
=  1.0  at  the  upper  surface  GHIJKL.  The  grid 
spacing  adjacent  to  the  surface  satisfied  Ayu */vw 
<  1.0.  A  uniform  z-grid  spacing  Az/8..  =  0.6 
was  employed.  The  computation  employing  the 
Baldwin-Lomax  model  utilized  a  grid  of  176,640 
points.  The  x-grid  spacing  was  concentrated  near 
the  comer,  with  a  minimum  Ax/ 8..  »  0.067. 
The  maximum  Ax/5..  *  0.2  and  0335  at  the 
upstream  and  downstream  boundaries,  respec¬ 
tively.  The  y-grid  spacing  was  concentrated  near 
the  surface  with  the  minimum  Ay/6..  * 

3.4x10'^,  and  the  maximum  Ay/6..  =  0.467  at 
the  upper  surface  GHIJKL.  The  grid  spacing 
adjacent  to  the  surface  satisfied  Ayu */vw  <  2.0. 
The  z-grid  spacing  was  concentrated  near  the 
symmetry  plane  AFEKLG,  with  a  minimum 
Az/5..  =  7.2x10"^,  and  a  maximum  Az/6*,  ■ 
0.564  near  the  boundary  BCDJIH. 

Comparison  with  Experiment 

R*&~  = 

Experimental  profiles  of  pitot  pressure  and 
yaw  angle  were  obtained  by  Settles  and  McKen¬ 
zie'^  at  fifteen  positions  for  the  (a,X)  =  ( 24,40) 
swept  compression  comer  at  Re&*  =  2.6xl(H, 
The  locations  of  the  measurements  are  shown  in 
Fig.  4.  The  profiles  were  taken  at  a  fixed  span- 
wise  position  z  =  9.686..,  and  the  stream  wise 
positions  (x-xcomeT)/ S„  =  -8.48,  -4.85,  -4.24, 
-3.64,  -2.73,  -0.61,  0.  0.61,  1.52,  3.03,  4.24, 
7.27,  9.09,  10.9,  and  14.6,  where  x^j^j  =  z 
tanX  is  the  location  of  the  comer  at  a  given  z. 
These  locations  are  designated  as  Positions  1,  2, 
....  15.  The  profiles  for  Positions  1  through  9  and 


15  were  obtained  along  the  vertical  axis.  The 
profiles  for  Positions  10  through  14  were  taken 
along  a  line  tilted  upstream  at  an  angle  y  =  16 
deg. 

A  set  of  four  pitot  pressure  and  yaw  angle 
profiles  have  been  selected  for  presentation.  The 
position  locations  were  chosen  based  upon  two 
criteria,  namely,  i)  the  comparison  between  the 
computed  and  experimental  profiles  is  indicative 
of  the  full  set  of  fifteen  profiles,  and  ii)  the  posi¬ 
tions  were  located  near  identifiable  features  of 
the  flow  (i.e.,  the  line  of  coalescence,  the  comer 
line,  the  line  of  divergence,  and  a  position  far 
downstream).  The  pitot  profiles  at  Position  5, 
located  approximately  at  the  line  of  coalescence, 
are  displayed  in  Fig.  5.  The  horizontal  axis  is  the 
pitot  pressure  Pp  normalized  by  the  upstream 
freestream  pitot  pressure  p_*, .  The  vertical  axis 
is  the  distance  measured  from  the  surface  (see 
above),  normalized  by  the  upstream  boundary 
layer  thickness  6M .  The  experimental  profile 
displays  an  approximate  40%  overshoot  in  Pp  at 
the  edge  of  the  boundary  layer.  This  overshoot 
is  associated  with  the  compression  region 
upstream  of  the  interaction.^"^  It  is  evident  to 
a  lesser  extent  in  the  calculations  employing  the 
Cebeci-Smith  and  Jones-Launder  models,  and  is 
absent  from  Baldwin-Lomax  computation.  The 
underestimate  of  the  pitot  pressure  overshoot  is 
attributable  to  the  underprediction  of  the 
upstream  influence  of  the  interaction  by  all  three 
models.  This  underprediction  is  also  apparent  in 
the  computed  surface  pressure  profiles  (see 
below). 

The  pitot  profiles  at  Position  7,  located  at 
the  comer  line,  are  shown  in  Fig.  6.  Good 
agreement  is  observed  between  the  computed 
and  measured  profiles,  with  discrepancies  limited 
to  the  vicinity  of  the  shock  wave  at  y  =  236*, 
where  the  computed  profiles  underpredict  the 
height  of  the  shock  wave.  The  underprediction 
is  another  manifestation  of  the  underestimate  of 
the  upstream  influence  of  the  interaction  (i.e., 
the  computed  shock  structure  is  situated  down¬ 
stream  of  the  experimental  shock  structure). 
The  pitot  profiles  at  Position  10,  located  approxi¬ 
mately  at  the  line  of  divergence,  are  displayed  in 
Fig.  7.  The  three  computed  profiles  are  very 
similar,  and  in  good  agreement  with  experiment. 
The  pitot  profiles  at  Position  14,  located  far 
downstream,  are  shown  in  Fig.  8.  Good  agree¬ 
ment  between  the  calculated  and  measured 


profiles  is  again  observed. 

The  yaw  angle  profiles  at  the  same  four 
Positions  are  displayed  in  Figs.  9  to  12.  The  yaw 
angle  f$  is  defined  as  (5=  tan‘*(w/(u  cosy  +  v 
siny)),  where  (u.v.w)  are  the  Cartesian  velocity 
components  in  the  (x,y,z)  coordinate  system, 
and  y  is  the  angle  of  tilt  of  the  measured  profile 
(see  above).  Near  the  line  of  coalescence  (Fig. 
9),  the  yaw  angle  is  observed  to  be  small  except 
in  the  immediate  vicinity  of  the  surface.  The 
computed  profiles  display  significant  differences 
near  the  wall.  The  proximity  of  this  position  to 
the  line  of  coalescence  accentuates  the 
differences  between  the  three  models,  since  the 
near-surface  yaw  angle  at  a  fixed  height  is 
observed  to  increase  rapidly  as  the  line  of  coales¬ 
cence  is  approached  from  upstream.  The  yaw 
angle  profiles  at  the  corner  line  (Fig.  10)  show 
close  agreement  between  the  three  models.  The 
experimental  data  is  limited  to  y  >  8„  due  to 
;he  size  of  the  boundary  layer  and  the  probe 
geometry.  Near  the  line  of  attachment  (Fig.  11), 
the  computed  profiles  are  nearly  identical,  and 
display  generally  good  agreement  with  experi¬ 
ment.  Finally,  the  computed  and  measured  yaw 
angle  profiles  far  downstream  (Fig.  12)  show 
close  agreement. 

The  computed  and  measured  surface  pres¬ 
sure  profiles  at  z  =  10.15..  are  shown  in  Fig.  13. 
There  is  general  agreement  between  the  calcu¬ 
lated  and  experimental  profiles,  although  the 
data  is  limited  to  points  upstream  of  the  corner 
line.  The  computations  using  the  Baldwin- 
Lomax  model  modestly  underpredict  the 
upstream  influence  (Le.,  the  location  of  the  ini¬ 
tial  pressure  rise).  Additional  data  at  t  > 
14.65..  exhibit  a  similar  behavior.  The  calcula¬ 
tions  utilizing  the  Cebeci-Smith  and  Jones- 
Launder  models  show  closer  agreement  with 
experiment.  These  computations  utilized  a 
coarser  grid,  however,  and  further  calculations 
are  needed  to  determine  the  effect  of  grid 
refinement. 

Reg^  =  8.1x10s 

Experimental  profiles  of  pitot  pressure  and 
yaw  angle  were  obtained  by  Settles  and  McKen¬ 
zie  at  fifteen  stations  for  the  (a, A.)  =  (24,40) 
swept  compression  comer  at  Reg^  =  8.1x10  . 
The  locations  of  the  measurements  are  shown  in 
Fig.  4.  The  profiles  were  taken  at  a  fixed  span- 


wise  position  z  =  75..,  and  streamwise  positions 

(*-*comer>/S~  =  ’4-  *3-2-  - 2 •».  '2-4.  -2.  -M. 
-0.8,  -0.2,  1.4,  2.6,  3.  3.6,  4.2,  4.8,  and  6.  Exper¬ 
imental  profiles  of  pitot  pressure,  yaw  angle  and 
static  pressure  were  obtained  by  Ketchum  and 
Bogdonoff  for  the  same  configuration  at  z  = 
75_,  and  streamwise  positions  (x-x-__.r)/5„ 
=  -4.0,  -2.8,  -1.4,  -0.8,  0.4,  1.6,  2.4,  3.6,  4.8,  6.0 
and  8.4.  The  resolution  of  these  profiles  (i.e., 
the  number  of  experimental  data  points  within 
the  boundary  layer)  was  greater  than  the  previ¬ 
ous  measurements  of  Settles  and  McKenzie. 
Seven  of  the  positions  of  Ketchum  and  Bog¬ 
donoff  coincide  with  the  Settles  and  McKenzie 
data.  The  two  sets  of  data  display  close  agree¬ 
ment  at  these  positions,  thereby  verifying  the 
repeatability  of  the  experiments  which  were  per¬ 
formed  several  years  apart.  The  two  data  sets 
have  been  merged  into  a  single  data  set  of  nine¬ 
teen  positions  at  (*'xcomerW®«  *  _4,  "3-2, 

-2.8,  -2.4,  -2,  -1.4,  -0.8,  -0.2,  0.4,  1.4,  1.6,  2.4, 
2.6,  3.0,  3.6,  4.2,  4.8,  6.0,  and  8.4.  The  profiles 
for  Positions  1  through  9,  11  and  12  were 
obtained  along  the  vertical  axis.  The  profiles  for 
Positions  10  and  13  through  17  were  obtained 
along  a  line  tilted  upstream  at  an  angle  f  «  16 
deg.  The  profiles  for  Positions  18  and  19  were 
taken  at  y  =  24  deg.  The  experimental  line  of 
coalescence13  is  observed  at  x-x^^g,.  =  -2.58„. 
for  Z  *  75„ ,  located  between  Positions  3  and  4. 
The  experimental  line  of  divergence  is  found  at 
x-Xcomer  =  2.65,,  for  z  =  75„,  located  at 

Position  13. 

A  set  of  four  pitot  pressure  and  yaw  angle 
positions  have  been  selected  for  presentation 
using  the  same  criteria  employed  for  Reg„  = 
2.6 x  10  ’,  In  Fig.  14,  computed  and  experimental 
pitot  pressure  profiles  are  displayed  at  Position  4. 
located  slightly  downstream  of  the  experimental 
line  of  coalescence.  The  profiles  are  in  close 
agreement,  and  display  negligible  effect  of  the 
interaction.  The  computed  and  measured 
profiles  at  Position  9,  located  0.45—  downstream 
of  the  comer,  are  displayed  in  Fig.  15.  The  cal¬ 
culated  profiles  agree  closely  with  the  experimen¬ 
tal  data  within  the  boundary  layer,  with  the 
Baldwin-Lomax  model  providing  a  slightly 
improved  prediction.  Outside  the  boundary 
layer,  the  computated  profiles  display  a  diffused 
profile  in  the  vicinity  of  the  shock  wave  due  to 
the  shock-capturing  nature  of  the  numerical  algo¬ 
rithms.  The  computed  and  experimental  pitot 
pressure  profiles  at  Position  13,  located  at  the 


Flowfield  Structure 


line  of  divergence,  are  shown  in  Fig.  16.  The 
profiles  display  excellent  agreement,  and  show 
negligible  variation  due  to  the  turbulence  model. 
Finally,  the  calculated  and  measured  profiles  at 
Position  18,  located  downstream  of  the  line  of 
divergence,  are  shown  in  Fig.  17.  Close  agree¬ 
ment  is  again  observed  between  the  computed 
and  measured  profiles. 

The  yaw  angle  profiles  at  the  same  four 
positions  are  displayed  in  Figs.  18  to  21.  The 
computed  and  measured  yaw  angle  at  Position  4, 
located  slightly  downstream  of  the  line  of  coales¬ 
cence,  display  negligible  yaw  angle  except  within 
the  immediate  vicinity  of  the  surface.  There  is 
considerable  discrepancy  between  the  calculated 
and  measured  surface  yaw  angles  at  this  position 
due  to  the  differences  in  the  calculated  upstream 
influence  (see  below).  The  calculated  and  exper¬ 
imental  yaw  angles  at  Position  9,  located  0.45„ 
downstream  of  the  comer,  are  shown  in  Fig.  19. 
The  computed  yaw  angle  profiles  are  virtually 
identical  indicating  insensitivity  to  the  tur¬ 
bulence  model  employed.  The  computed  profiles 
overpredict  the  experimental  yaw  angle  profile  of 
Ketchum  and  Bogdonoff  in  the  outer  portion  of 
the  boundary  layer.  The  computed  and  meas¬ 
ured  yaw  angle  profiles  at  Position  13,  located  at 
the  line  of  divergence,  are  displayed  in  Fig.  20. 
The  calculated  profiles  are  seen  to  be  virtually 
identical,  and  in  close  agreement  with  experi¬ 
ment.  Finally,  the  yaw  angle  profiles  at  Position 
18.  located  downstream  of  the  line  of  diver¬ 
gence,  are  presented  in  Fig.  21.  Excellent  agree¬ 
ment  is  observed  between  the  calculated  profiles 
and  experiment. 

The  computed  and  measured  surface  pres¬ 
sure  profiles  at  z  *  105„  are  shown  in  Fig.  22. 
The  calculated  profile  using  the  Baldwin-Lom ax 
model  underpredicts  the  upstream  influence  (i.e., 
the  location  of  the  initial  pressure  rise).  This 
behavior  has  been  observed  previously  in  a 
variety  of  3-D  shock  wave-turbulent  boundary 
layer  computations  using  the  Baldwin-Lomax 
model.11'14  The  plateau  pressure  is  accurately 
predicted.  The  calculated  profiles  using  the 
Cebeci-Smith  and  Jones-Launder  models  provide 
a  more  accurate  prediction  of  the  upstream 
influence,  but  significantly  overpredict  the  pla¬ 
teau  pressure.  Further  investigation  of  the  effect 
of  grid  refinement  for  the  calculations  using  the 
Cebeci-Smith  and  Jones-Launder  models  is 
required. 


The  general  agreement  between  computed 
and  experimental  pitot  pressure  and  yaw  angie 
profiles  for  the  (a,X)  *  (24,40)  deg 

configuration  (see  above)  and  the  (a.X)  = 
(24,60)  deg  configuration14  permits  the  develop¬ 
ment  of  flowfield  models  based  on  the  computed 
results.  The  calculated  profiles  of  pitot  pressure 
and  yaw  angle  have  been  observed  to  be  insensi¬ 
tive  to  the  turbulence  model  employed,  except 
within  the  immediate  vicinity  of  the  surface 
(e.g.,  typically  the  lower  10%  of  the  boundary 
layer).  It  suffices,  therefore,  to  examine  the 
flowfield  structure  using  the  computed  flowfields 
of  a  single  turbulence  model  (e.g.,  Baldwin- 
Lomax),  with  the  understanding  that  some  quan¬ 
titative  differences  will  appear  near  the  surface. 
The  present  paper,  therefore,  focuses  on  a 
detailed  quantitative  study  of  the  flowfield  struc¬ 
ture  for  the  (a ,X)  =  (24,40)  and  (24,60)  deg 
configurations  using  the  computed  flowfields  with 
the  Baldwin-Lomax  model. 

The  computed  surface  skin  friction  lines 
for  the  (a,X)  *  (24.40)  configuration  are 

displayed  in  Figs.  23  and  24  for  Reg^  * 
2.6x10^  and  8.1x10^.  respectively.  The  line  of 
coalescence,  formed  by  the  convergence  of  the 
upstream  skin  friction  lines,  is  clearly  evident. 
Results  for  the  (a.X)  «  (24,60)  configuration14 
indicated  that  angle  of  the  computed  line  of 
coalescence,  measured  relative  to  the  span  wise 
direction  z,  was  approximately  10%  greater  than 
the  experimental  value  based  upon  kerosene 
lampblack  surface  visualization.  The  line  of 
divergence,  situated  on  the  compression  surface, 
is  also  shown. 

The  flowfield  structures  for  the  (a,X)  <= 
(24,40)  deg  and  (24,60)  deg  configurations  are 
similar.  A  general  mean  flowfield  model, 
developed  on  the  basis  of  particle  pathlines14,  is 
presented  in  Fig.  25.  It  is  emphasized  that  this  is 
a  model  based  upon  the  mean  flowfield.  whereas 
in  the  actual  experiment  the  flowfield  is  tur¬ 
bulent  and  therefore  unsteady.  The  principal 
element  of  the  flow  is  a  large  vortical  structure. 
The  line  of  coalescence  defines  the  origin  of  a 
three  dimensional  surface  of  separation,  denoted 
as  Surface  No.  1.  The  surface  spirals  into  the 
core  of  the  vortical  structure.  It  should  be 
noted,  however,  that  the  streamlines  comprising 
Surface  No.  1  exhibit  a  significant  yaw  angle,  and 
consequently  are  skewed  strongly  in  the  span  wise 
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direction.  Approaching  the  line  of  coalescence, 
the  streamlines  close  to  the  surface  rise 
moderately  and  move  strongly  in  the  spanwise 
direction.  At  the  computed  line  of  coalescence 
for  z  =  75^.  for  example,  the  maximum  angle 
of  pitch  (defined  as  tan'*(v/V  u^+  w^)  is  approx¬ 
imately  4  deg  and  occurs  at  y  =  0.125^.;  the 
maximum  yaw  angle  is  60  deg  and  occurs  at  the 
surface.  The  line  of  divergence  defines  the  inter¬ 
section  of  a  second  surface,  denoted  as  Surface 
No.  2,  with  the  swept  compression  ramp.  This 
second  surface  extends  upstream  into  the  undis¬ 
turbed  flow.  The  fluid  contained  between  the 
wall  and  Surface  No.  2  becomes  the  vortical 
structure,  which  is  approximately  aligned  with 
the  corner.  The  fluid  above  the  second  surface 
passes  over  the  vortical  structure  and  continues 
up  the  compression  ramp. 

A  quantitative  measure  of  the  Surface  No. 

1  is  its  angle  of  intersection  with  the  wall,  meas¬ 
ured  in  a  plane  orthogonal  to  the  line  of  coales¬ 
cence.  The  calculated  angle  is  10  t.  2  deg  for  the 
(0.7)  *  (24,40)  deg  configuration  at  Rej^  * 
2.6x10^  and  8.1x10^.  The  same  value  is 
obtained  for  the  (24,60)  deg  configuration  at 
Rej^  m  1.4x10^  and  9x10^.  The  angle  of  inter¬ 
section  is  not,  however,  the  angle  of  pitch  of  the 
streamlines  comprising  Surface  No.  1.  Those 
streamlines  are  moving  spanwise  in  a  direction 
approximately  aligned  with  the  line  of  coales¬ 
cence,  and  therefore  their  pitch  is  less  than  the 
angle  of  intersection  of  Surface  No.  1  with  the 
wall. 

An  important  feature  of  Surface  No.  2  is 
its  height  above  the  wall,  measured  at  a  fixed 
stream  wise  position  immediately  upstream  of  the 
interaction.  The  calculated  height  of  Surface  No. 

2  is  displayed  in  Fig.  26  for  (a,y)  ■  (24,40)  deg 
at  ^e6«.  *  2.6xl03  and  8.1  xlO3,  and  for  (a,X) 
»  (24,60)  deg  configuration  at  Reg^  =  1.4x10* 
and  9x10^.  The  height  of  Surface  No.  2  at  a 
given  spanwise  location  z  is  determined  by  exa¬ 
mining  streamlines  originating  at  various  values 
of  y,  and  observing  whether  the  particular 
streamline  moves  up  along  the  compression  sur¬ 
face  or  spirals  into  the  vortical  structure.  The 
height  is  defined  as  the  limit  where  the  stream¬ 
lines  approach  the  line  of  divergence  and  move 
parallel  to  it.  As  the  value  of  z  increases,  greater 
uncertainty  arises  in  determining  whether  a 
streamline  is  moving  up  the  ramp  or  mto  the 
vortex.  This  increasing  uncertainty  is  associated 


with  the  limited  spanwise  extent  of  the  computa¬ 
tional  domain  to  view  the  particular  streamline. 
Consequently,  there  are  two  lines  shown  for  each 
value  of  (a,X)  and  Reg^ .  These  lines  represent 
the  upper  and  lower  estimates  of  the  height  of 
Surface  No.  2. 

There  are  two  important  features  regarding 
the  height  of  Surface  No.  2.  First,  the  height  for 
the  (a,X)  =  (24,40)  deg  configuration  appears  to 
asymptote  at  y  *  0.155.,,  while  for  the  (a.X.) 
=  (24,60)  deg  configuration  the  height  continues 
to  increase  in  an  approximate  linear  manner  well 
beyond  y  «  5„.  The  limited  spanwise  extent  of 
the  computations  precludes  a  definitive  statement 
regarding  the  asymptotic  behavior  of  the  charac¬ 
teristics  of  the  computed  flowfields  including  the 
height  of  Surface  No.  2  and  the  nature  of  the 
line  of  coalescence.  Settles  and  Teng^  classified 
both  the  (O.X)  -  (24.40)  and  (24,60) 

configurations  as  "conical"  in  the  sense  that  the 
line  of  coalescence  was  determined  to  asymptote 
a  straight  line  at  a  fixed  angle  relative  to  the 
comer  line,  whereas  Wang  and  Bogdonoff^  sug¬ 
gested  that  the  farfield  behavior  would  always 
approach  cylindrical  flow.  The  relationship 
between  the  asymptotic  behavior  of  the  line  of 
coalescence  and  the  height  of  Surface  No.  2  is 
not  understood  at  present,  and  is  the  topic  of 
further  investigation.  Second,  the  results  for  the 
two  different  Reynolds  numbers  at  each  (a.X.) 
indicate  that  the  height  appears  to  scale  with  the 
upstream  boundary  layer  thickness.  McClure 
and  Dolling^3  suggested  that  the  flowfield  struc¬ 
ture  of  the  3-D  sharp  fin  at  o  =  10  deg,  Mach  3 
and  Rej^  *  3.2xlCH  and  7.9x10^  scaled  accord¬ 
ing  to  (y/50)Re5o1/3  vs.  (i/S^Re^,"1/3 
based  upon  surface  measurements  of  the  3-D 
swept  compression  comer  at  similar  flow  condi¬ 
tions  by  Settles  and  Bogdonoff^,  where  50  is 
the  local  incoming  boundary  layer  thickness. 
The  present  computations,  although  limited  in 
spanwise  extent,  do  not  support  the  flowfield 
scaling  model  of  McClure  and  Dolling. 

Further  understanding  of  the  flowfield 
structure  may  be  obtained  through  examination 
of  the  governing  equations,  specifically,  the 
evaluation  of  the  various  terms  in  the  governing 
equations  along  a  streamline.  The  mean  kinetic 
energy  per  unit  mass  ("mechanical  energy")  is 
examined  for  simplicity.  The  equation  is 

D(Vto|Uj)/Dt  =  -(uj/p)3p/3xj  ♦  (Uj/pldtjj/dxj 
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where  ViujUj  is  the  mean  kinetic  energy  per  unit 
mass,  D/Dt  *  3/3t  +  Uj9/3xj,  (uj.u2.u3)  * 
(u,v,w).  (xj.x2.x3)  “  (x.y.z).  p  is  the  static 

pressure,  p  is  the  density,  Tjj  is  the  total  (viscous 
plus  turbulent)  stress  tensor,  and  the  Einstein 
summation  convention  is  employed.  There  are 
two  contributions  to  the  rate  of  change  of  the 
mean  kinetic  energy  per  unit  mass.  The  first 
term.  -(uj/p)3p/3xj,  is  an  inviscid  effect  associ¬ 
ated  with  compression  or  dilation.  The  second 
term,  (Uj/p)9tjj/3xj.  represents  the  effect  of  the 
viscous  and  turbulent  stresses.  A  series  of  flow 
streamlines  are  generated  by  integrating  the 
motion  of  a  fluid  particle,  released  at  a  selected 
upstream  position,  through  the  computed 
flowfield.  The  mean  kinetic  energy  equation  is 
integrated  simultaneously  along  the  streamline. 


V4ujuj  =  ‘Aujujlo  + 

j  l-(ui/p)dp/dxi  +  (Uj/pldtij/axjjds/v 

where  V4ujUj|0  is  the  initial  mean  kinetic  energy 
at  the  point  of  release  of  the  fluid  particle,  ds  is 
the  infinitesimal  arclength  along  the  streamline 
and  V  *  V  ujuj.  The  separate  contributions  of 
the  integrands  -(uj/p)3p/3xj  and  (uj/p)dxjydxj 
provide  a  quantitative  measure  of  the  inviscia 
and  viscous  contributions  to  the  rate  of  change 
of  the  mean  kinetic  energy. 

The  contributions  to  the  rate  of  change  of 
the  mean  kinetic  energy  are  evaluated  for  the 
(a.X)  =  (24,40)  deg  configuration  at  Re$_  = 
8.1x10^  from  the  calculated  flowfield  obtained 
using  the  Baldwin-Lomax  model  A  series  of  ten 
streamlines,  released  within  the  upstream  undis¬ 
turbed  boundary  layer  at  y  »  0.028„ ,  are 

displayed  in  Fig.  27.  The  entrainment  of  the 
streamlines  into  the  vortical  structure  is  clearly 
evident.  The  mean  kinetic  energy  on  the  fifth 
streamline  (counted  from  the  plane  of  sym¬ 
metry)  is  presented  in  Fig.  28  as  a  function  of 
arclength  along  the  streamline.  The  figure 
displays  both  the  integrated  kinetic  energy, 
obtained  from  the  path  integral  of  the  mean 
kinetic  energy  equation  (see  above),  and  the 
interpolated  kinetic  energy,  obtained  from  inter¬ 
polation  of  the  computed  flowfield.  These  two 
results  are  not  identical  due  to  the  truncation 
errors  of  the  numerical  computation  of  both  the 
flowfield  and  the  path  integral.  The  results,  how¬ 
ever,  are  seen  to  be  in  close  agreement.  The 
kinetic  energy  decreases  abruptly  at  the  start  of 


the  interaction,  reaches  a  minimum  value  shortly 
after  entering  the  vortical  structure,  and  eventu¬ 
ally  asymptotes  a  nearly  constant  value.  In  Fig. 
29,  two  profiles  are  displayed  for  this  streamline. 
The  first  profile,  denoted  the  "Pressure  and 
Viscous  Contribution",  is  the  integrated  mean 
kinetic  energy  displayed  previously  in  Fig.  28, 
le,  the  path  integral  of  the  mean  kinetic  energy 
equation  obtained  by  including  both  the  pressure 
and  viscous  integrands ,  i.e„  -(uj/p)3p/3xj  and 
(Uj/p)dTjj/dxj.  The  second  profile,  denoted  the 
Treasure  Contribution",  is  the  path  integral  of 
the  mean  kinetic  energy  equation  obtained  by 
including  only  the  pressure  integrand,  i.e., 
-(Uj/p)3p/dxj.  The  two  profiles  are  observed  to 
be  significantly  different,  and  it  is  therefore  con¬ 
cluded  that  both  inviscid  and  viscous  effects  are 
important  for  this  particular  streamline.  Similar 
results  are  observed  for  the  other  nine  stream¬ 
lines  emanating  at  y*  0.025... 

In  Fig.  30,  a  series  of  ten  streamlines, 
released  at  y  *  0.055..,  are  displayed.  The 
streamlines  roll  up  into  the  vortical  structure  as 
expected  from  Fig.  26.  The  interpolated  and 
integrated  kinetic  energy  for  the  fifth  streamline, 
displayed  in  Fig.  31,  are  in  close  agreement.  The 
profiles  in  Fig.  32  display  the  mean  kinetic 
energy  of  the  fifth  streamline  obtained  from  the 
path  integral  of  the  pressure  and  viscous  contri¬ 
butions,  and  from  the  path  integral  of  the  pres¬ 
sure  contribution  alone.  The  integrated  mean 
kinetic  energy,  obtained  from  the  path  integral  of  the 
pressure  contribution  alone,  provides  a  reasonably 
accurate  measure  of  the  kinetic  energy.  This 
implies,  therefore,  that  the  turbulent  stresses 
have  a  minimal  effect  on  the  evolution  of  these 
streamlines,  despite  their  origin  deep  within  the 
boundary  layer  at  y  *  0.05  6..  and  their  entrain¬ 
ment  within  the  vortical  structure.  Similar  results 
are  observed  for  the  other  nine  streamlines 
emanating  at  y  *  0.055„. 

In  Fig.  33,  the  final  series  of  ten  stream¬ 
lines,  released  at  y  *  0.55.. ,  are  shown.  These 
streamlines  move  over  the  vortical  structure,  and 
continue  up  along  the  ramp  (Fig.  26).  The 
integrated  and  interpolated  mean  kinetic  energy 
for  the  fifth  streamline,  displayed  in  Fig.  34, 
show  a  discrepancy  of  approximately  18%. 
Further  investigations  are  in  progress  to  examine 
this  effect.  The  profiles  in  Fig.  34  display  the 
mean  kinetic  energy  of  the  fifth  streamline 
obtained  from  the  path  integral  of  the  pressure 
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and  viscous  contributions,  and  the  path  intergral 
of  the  pressure  contribution  alone.  It  is  again 
observed  that  the  mean  kinetic  energy  obtained 
from  the  path  integral  of  the  pressure  contribu¬ 
tion  alone  is  a  reasonable  prediction,  varying 
from  the  full  interpolated  kinetic  energy  by  only 
15%.  Similar  results  are  observed  for  the  other 
nine  streamlines  emanating  at  y  *  0.56...  It  is 
again  evident  that  turbulent  stresses  play  a  minor 
role  in  the  dynamics  of  the  fluid  motion  for 
these  streamlines. 

The  examination  of  the  mean  kinetic 
energy  for  the  (a.X)  ■  (24,40)  deg  configuration 
at  Re  5..  *  8.1x10^  provides  a  quantitative  indi¬ 
cation  of  the  importance  of  inviscid  and  viscous 
effects.  Except  in  the  immediate  vicinity  of  the 
surface,  the  flowfield  is  observed  to  be  rotational 
and  inviscid,  as  conjectured  previously  for  the 
3-D  sharp  fin.^  Further  investigations  are  in 
progress  for  other  swept  compression  comer 
configurations. 

Conclusions 

A  collaborative  experimental  and  theoreti¬ 
cal  research  program  has  focused  on  the  3-D 
shock  wave-turbulent  boundary  layer  interaction 
generated  by  a  swept  compression  comer, 
characterized  by  the  compression  angle  a  and  the 
sweep  angle  X.  In  the  present  effort,  computa¬ 
tions  have  been  performed  for  (a,X)  -  (24,40) 
deg  at  Rej^  *  2.6x10^  and  8.1x10^  at  Mach  3 
using  the  Baldwin-Lomax  turbulence  model. 
These  results  are  compared  with  the  experimen¬ 
tal  data  of  Settles  and  McKenzie,  Ketchum  and 
Bogdonoff,  and  previous  computations  by  Horst- 
man  using  the  Cebeci-Smith  and  Jones-Launder 
turbulence  models.  The  calculated  and  experi¬ 
mental  profiles  of  pitot  pressure  and  yaw  angle 
are  in  good  agreement.  The  computed  profiles 
are  remarkably  insensitive  to  the  turbulence 
model  employed  except  in  a  narrow  region 
comprising  the  lower  10%  of  the  boundary  layer 
where  differences  of  10%-20%  in  yaw  angle  are 
observed.  The  calculated  surface  pressure  using 
the  Baldwin-Lomax  underpredicts  the  location  of 
upstream  influence,  although  providing  reason¬ 
able  agreement  elsewhere.  The  computed  sur¬ 
face  pressure  using  the  Cebeci-Smith  and  Jones- 
Launder  models  provides  a  better  prediction  of 
the  upstream  influence  point,  but  overpredict  the 
pressure  in  the  vicinity  of  the  comer.  Further 
calculations  using  the  Cebeci-Smith  and  Jones- 
Launder  models  are  required  to  ascertain  the 


The  dominant  feature  of  the  (a,X)  = 
(24,40)  deg  flowfield  is  a  large  vortical  structure, 
approximately  aligned  with  the  comer,  similar  to 
the  flow  structure  for  the  (a.X)  =  (24,60)  deg 
configuration.  A  surface  of  streamlines 
emanates  from  the  line  of  coalescence,  and 
forms  the  core  of  the  vortical  structure.  The 
surface  intersects  the  wall  at  an  angle  of  10-t2 
deg  for  both  configurations.  The  angle  of  inter¬ 
section  is  not,  however,  the  angle  of  pitch  of  the 
streamlines  comprising  Surface  No.  1.  Those 
streamlines  are  moving  span  wise  in  a  direction 
approximately  aligned  with  the  line  of  coales¬ 
cence,  and  their  pitch  is  therefore  less  than  the 
angle  of  intersection  of  Surface  No.  1  with  the 
wall.  A  second  surface  originates  upstream 
within  the  boundary  layer  and  interseas  the 
compression  ramp  at  the  line  of  divergence.  The 
fluid  contained  between  the  wall  and  the  second 
surface  is  entrained  into  the  vortical  structure. 
The  height  of  the  surface,  measured  immediately 
upstream  of  the  interaction,  is  strongly  depen¬ 
dent  on  the  sweepback  angle  X  and  scales  with 
the  upstream  boundary  6..  . 

The  mean  kinetic  energy  equation  is  exam¬ 
ined  for  a  series  of  streamlines  emanating  within 
the  upstream  boundary  layer  for  the  (a.X)  * 
(24,40)  deg  configuration  at  Rej^  «  8.1xl05. 
Three  initial  heights  were  selected,  namely, 
y/8„  *  0.02,  0.05  and  0.5.  The  change  of 

kinetic  energy  along  streamlines  emanating  from 
y  «  0.025„  was  observed  to  be  influenced 

strongly  by  both  inviscid  and  viscous  effeas. 
However,  for  the  streamlines  originating  at  y/8— 
=  0.05  and  0.5,  the  rate  of  change  of  kinetic 
energy  was  attributable  almost  entirely  to  inviscid 
effects.  The  streamlines  emanating  at  y  * 
0.05S„  were  entrained  within  the  vortical  struc¬ 
ture,  while  the  streamlines  originating  at  y  = 
0.58..  moved  over  the  vortical  structure  and  up 
along  the  compression  surface.  This  represents 
the  first  direct  indication  that  the  vortical  struc¬ 
ture  is  predominantly  rotational  and  inviscid. 
except  within  a  very  narrow  region  of  the  boun¬ 
dary  layer  adjacent  to  the  wall. 
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ABSTRACT 

A  numerical  investigation  is  focussed  on  the  effect  of 
bleed  on  the  three-dimensional  shock  wave  -  turbulent 
boundary  layer  interaction  generated  by  a  3-D  sharp  fin 
with  20*  wedge  angle  at  a  Reynolds  number  9x10s  based 
on  the  upstream  boundary  layer  thickness.  The  Reynolds- 
averaged  3-D  compressible  Navier-Stokes  equations  in 
mass  averaged  variables  are  solved  with  the  eddy 
viscosity  as  prescribed  by  the  Baldwin- Lomax  algebraic 
model.  Bleed  is  applied  in  the  region  between  the  line  of 
upstream  influence  and  the  theoretical  inviscid  shock  line 
in  an  attempt  to  reduce  or  eliminate  altogether  the 
undesirable  separation  observed  in  such  flows.  Three 
magnitudes  of  mass  flux  (1%,  2.5%  and  5%  of  free 
stream)  are  computed.  The  computed  flow  fields  are 
compared  with  experimental  and  theoretical  results 
obtained  in  the  absence  of  bleed.  The  effect  of  the 
described  bleed  schedule  on  the  flow  field  is  modest  and 
relatively  local.  Various  physical  variables  such  as  the 
pitot  pressure  and  yaw  angles  are  moderately  influenced 
only  at  dose  proximity  of  the  surface.  In  the  interaction 
region,  the  previously  observed  overshoot  in  pitot  pressure 
profile  is  reduced.  Skin  friction  values  increase 
significantly  only  in  the  region  of  bleed.  The  line  of 
coalescence  is  seen  to  align  with  the  theoretical  inviscid 
shock  line  though  this  may  be  due  to  the  particular 
region  of  bleed  under  consideration.  A  detailed  study  of 
particle  pathlines  generated  from  the  computed  flow  field 
indicates  that  the  effect  of  bleed  is  essentially  to  ingest  a 
portion  of  the  boundary  layer.  The  vortical  structure 
observed  in  previous  research  is  not  significantly  affected. 
The  paper  also  describes  results  of  a  numerical 
investigation  of  a  weak  shock-shock  intersection  in  the 
presence  of  a  turbulent  boundary  layer.  Previous  research 
indicates  that  such  interactions  may  be  capable  of 
producing  a  given  pressure  rise  with  less  likelihood  of 
separation  than  an  equivalent  strength  single  shock 
interaction.  Numerical  results  of  flow  due  to  two  five 
degree  symmetrically  placed  wedges  on  a  flat  plate  are 
compared  with  experiment. 

X  Introduction 

3-D  shock  wave  -  turbulent  boundary  layer  interactions 
(denoted  "3-D  turbulent  interactions”)  occur  in  a  variety 
of  applications  such  as  aircraft  inlets  and  compressors, 
wing-body  junctures  and  control  surface  deflections.  In  the 
past  few  years  significant  progress  has  been  achieved  by 
theoretical  and  experimental  researchers  in  the 
understanding  of  3-D  turbulent  interactions  generated  by 
dimensionless  geometries  such  as  the  sharp  fin  (Fig.  1) 
which  consists  of  an  unswept  wedge  at  angle  of  attack  a 
mounted  on  a  flat  plate  on  which  develops  the  incoming 
turbulent  boundary  layer. 
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The  principal  parameters  in  such  flows  are  the  Mach 
number  Ma,  the  Reynolds  number,  Re^  ,  based  upon 

the  boundary  layer  thickness  l  at  the  streamwise  station 
corresponding  to  the  leading  edge  of  the  fin,  the  thermal 
boundary  conditions  and  the  fin  angle  a.  Experimental 
results  available  for  such  flows  may  be  classified  into 

surface  1,  2,  3,  4,  5,  6j  and  boundary  layer 

measurements  [7,  8,  9,  10,  llj.  These  experiments  cover 
a  range  of  the  above  mentioned  parameters  and  have 
provided  valuable  insight  into  the  structure  of  3-D 

turbulent  interactions. 

While  a  few  researchers  have  achieved  success  in  the 

analytical  treatment  of  shock  •  boundary  layer 

interactions  especially  for  weaker  shock 

strengths  12,  13,  14,  15J,  the  thrust  of  the  theoretical 
effort  is  mainly  numerical.  With  the  advent  of  high¬ 
speed  computers,  it  has  become  possible  to  numerically 
simulate  3-D  turbulent  interactions  '16,  17,  18,  191. 
These  computations  typically  employ  the  Reynolds 
averaged  Navier-Stokes  equations  coupled  with  an 
algebraic  20)  or  two-equation  (21)  turbulence  model  and 
have  demonstrated  good  agreement  with  experimental 
boundary  layer  and  surface  measurements  at  Mach  3  for 
Re4  =  2.5  x  10s  to  9.0  x  10s  and  wedge  angles  up  to 

OD 

20*  [10].  Results  obtained  computationally  for  the  swept- 
compression  corner  have  also  been  encouraging  [221. 

Recent  calculations  at  fin  angle  20*  by  Knight  et  ai  '  10) 
with  two  different  turbulence  models,  have  proven  to  be 
in  close  agreement  with  experiment  with  modest 
discrepancies  in  the  immediate  vicinity  of  the  flat  plate. 
Analysis  of  these  computations  has  led  to  the  conclusion 
that  the  principal  flow  feature  is  a  large  vortical 
structure  aligned  with  the  corner  in  agreement  with  the 
flow  field  models  of  Token  23]  and  Kubota  and 
Stollery  4i.  A  three-dimensional  surface  of  separation 
(Fig.  2)  emanates  from  the  line  of  coalescence 
(separation),  and  spirals  into  the  vortical  center.  A 
second  surface,  emanating  from  upstream,  intersects  the 
wail  at  the  line  of  divergence  (attachment),  and  defines 
the  extent  of  the  fluid  entrained  into  the  vortical 
structure. 

The  characteristics  of  the  above  described  flow  field  - 
specifically  the  large  vortical  structure  and  separated  flow 
-  are  not  ideal  in  applications  such  as  aircraft  inlets. 
Improvements  in  the  control  of  high'  speed  3-D  turbulent 
interactions  may  lead,  therefore,  to  improved  aircraft 
performance.  Several  possible  means  of  control  may  be 
identified: 

•  Bleed  (suction).  Bleed  has  been  traditionally 
employed  in  high  speed  aircraft  inlets  to 
prevent  boundary  layer  separation.  The  effect 
of  suction  is  to  remove  the  low  speed  fluid  in 
the  boundary  layer  before  it  separates  from  the 
surface. 
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•  Blowing:  Experimental  measurements  by 

Peake  8]  at  Mach  numbers  2  and  4  and 
Reynolds  number  2  x  10s  based  on  the 
undisturbed  boundary  layer  thickness  indicate 
that  it  is  possible  to  control  separation  by 
tangential  air  injection.  The  theory  behind  this 
approach  is  to  supply  additional  energy  to  the 
fluid  being  retarded  in  the  boundary  layer  by 
the  pressure  rise. 

•  Shock-Shock  Interaction:  Mae  and  Stalker  '151 
have  concluded  from  experimental  observations 
on  weak  shocks  that  intersecting  shocks  can 
produce  a  given  overall  pressure  rise  with  less 
likelihood  of  separation  than  an  equivalent 
strength  single  shock  interaction. 

•  Vorticitv  Distribution  Modification:  Sines  the 
principal  flow  feature  is  a  vortical  structure,  it 
may  be  possible  to  appiy  control  through  the 
introduction  of  additional  vorticity  in  the 
longitudinal  or  spanwise  direction. 

•  Modification  of  Geometry:  It  may  be  possible 
to  significantly  influence  the  flow  field  by 
minor  modifications  in  the  3-D  configuration 
itself.  The  effect  of  raising  the  fin  a  slight 
distance  from  the  flat  plate  could  prove 
interesting. 

The  objective  of  the  present  research  effort  is  to  study 
these  techniques  of  control  through  numerical  simulation. 
In  this  paper  the  first  of  these  rhoices,  namely  bleed 
(suction),  is  investigated.  The  numerical  simulation  of  a 
weak  shock-shock  interaction  in  the  presence  of  a 
turbulent  boundary  layer  is  also  described.  It  is 
recognised  that  there  is  a  scarcity  of  published  work 
describing  physical  experiments  incorporating  bleed  in  3-D 
turbulent  interactions  for  the  sharp  fin  configuration  and 
that  it  is  therefore  not  possible  to  make  direct 
comparison  of  the  results  to  be  presented  with 
experiment.  It  may  be  emphasised,  however,  that  the 
governing  Reynolds  averaged  3-D  compressible  Navier- 
Stokes  equations  and  the  numerical  algorithm  employed 
(the  hybrid  explicit-implicit  method  of  Knight  [17])  have 
been  applied  successfully  to  predict  turbulent  interactions 
for  the  3-D  sharp  fin  [10|  and  for  a  simulated  3-D  inlet 
with  bleed  [24].  In  addition,  the  algebraic  turbulent  eddy 
viscosity  model  of  Baldwin  and  Lomax  [201  has  been 
validated  for  a  variety  of  boundary  layer  flows  [251  with 
zero,  favorable  and  adverse  pressure  gradients  and  in  the 
presence  of  suction.  It  is  believed  therefore  that  the 
results  presented  serve  two  purposes,  namely,  to  improve 
the  current  state  of  understanding  of  3-D  turbulent 
interactions  and,  to  "reduce  the  required  experimental! 
test  matrix  to  the  smallest  number  of  configurations  26!" 
by  identifying  broad  trends  should  bleed  be  applied  to 
such  flows. 


2  Description  of  computations 


2.1  Governing  equations  and  numerical  model: 

The  governing  equations  are  the  full  mean  compressible 
Navier-Stokes  equations  using  mass  averaged  variables  and 
strong  conservation  form.  The  molecular  and  turbulent 
Prandtl  numbers  are  0.73  and  0.9,  respectively 


Turbulence  is  modeled  through  the  inclusion  of  the  two- 
layer  algebraic  turbulent  eddy  viscosity  model  of  Baldwin 
and  Lomax  20!  with  the  mixing  length  (f)  as  specified  by 
the  formula  of  Buieev  [271.  The  effect  of  bleed  is 
incorporated  in  the  Van  Driest  Damping  Factor  D  given 
by: 
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where  jr^;  is  the  magnitude  of  the  wall  shear  stress,  pm  is 
the  wall  dynamic  viscosity,  pm  is  the  fluid  density  at  the 
wall  and  N  is  the  bleed  correction  factor  of  Cebeci  28 
given  by 
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where  m  is  the  normal  mass  flux  at  the  wall,  with  m 
negative  for  bleed. 

A  3-D  coordinate  transformation  is  employed  as  in 
[10]  to  facilitate  the  application  of  the  numerical 
algorithm  and  the  boundary  conditions.  The  hybrid 
explicit-implicit  algorithm  of  Knight  [24]  is  employed  to 
solve  the  governing  equations.  The  implicit  algorithm 
employs  Keller’s  Box  Scheme  [29]  and  is  applied  to  the 
asymptotic  form  of  the  Navier-Stokes  equations  in  a  thin 
region  adjacent  to  the  solid  surfaces.  This  region,  denoted 
the  ’computational  sublayer’,  is  defined  by  z’+  <  SO, 
where  »’*  =  z’u./^,  with  u.  =  (rjpw)l/t,  rw  ="  wall 
shear  stress,  and  vm  is  the  kinematic  viscoeity  at  the 
wall.  Typically,  the  height  of  this  region  is  less  than  1% 
of  the  local  boundary  layer  thickness.  The  asymptotic 
form  of  the  conservation  of  momentum,  given  in  [17], 
represents  a  balance  between  convection  normal  to  the 
surface  (associated  with  bleed),  pressure  forces  and  viscous 
diffusion  normal  to  the  surface.  The  asymptotic  form  of 
the  conservation  of  energy  [17]  is  obtained  by  neglecting 
the  effects  of  stream  wise  convection  and  diffusion  and 
represents  a  balance  between  the  total  enthalpy 
transported  normal  to  the  surface  due  to  convection 
(associated  with  the  bleed)  and  viscous  diffusion.  The 
explicit  algorithm  of  MacCormack  301  is  applied  to  the 
full  Navier-Stokes  equations  on  grid  points  exterior  to  the 
computational  sublayer.  This  includes  nearly  all  of  the 
boundary  layer  and  the  external  inviscid  region.  The 
accuracy  of  the  computational  sublayer  approach  has  been 
verified  for  2-D  shock  wave-turbulent  boundary  layer 
interactions  in  the  presence  of  separation  [31],  and  is 
examined  in  Section  2.4  for  the  present  cases. 


3.3  Boundary  Conditions: 

The  boundary  conditions  are  (See  Fig.  1): 

•  Upstream  boundary  ABHG1:  The  flow  at  this 

boundary  is  assumed  to  be  2-D  and  is 

generated  with  a  separate  code  25!  such  that 
the  momentum  thickness  at  the  leading  edge 
of  the  fin  is  in  close  agreement  with 

experiments  of  Shapey  and  Bodgonoff  101.  It 
is  emphasized  that  these  experiments  do  not 
incorporate  suction.  The  various  flow 

parameters  are  indicated  in  Table  1. 

•  Plane  of  symmetry  AFLG  :  The  normal 
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component  of  the  velocity  end  the  normal 
derivetives  of  the  remaining  flow  quantities  ere 
set  to  sero. 

•  Richt  boundary  BCDJIH!:  Thij  boundary  is 
assumed  to  be  sufficiently  far  from  the  fin  to 
insure  that  the  boundary  layer  is  locally  two- 
dimensional  and  therefore  a  simple  gradient 
boundary  condition  d  ,dt— 0  is  employed. 

•  Downstream  boundary  EKJD  :  The 

conventional  sero  gradient  condition  is 
specified. 

•  Fin  surface  'LFKE!:  Since  this  is  a  solid 

surface,  the  velocity  and  the  normal  pressure 
gradient  are  taken  to  be  tero  and  a  fixed 
surface  temperature  (T_/T.  =  1.17)  is 

specified. 

•  Fiat  plate  lABCDEFI:  The  two  components  of 

velocity  along  the  flat  plate  are  specified  to  be 
sero.  Bleed  is  applied  in  the  triangular  area 

AYXA  shown  in  Fig.  3.  The  upstream 

boundary  AX  coincides  approximately  with  the 
line  of  upstream  influence,  defined  as  the  line 
where  the  incoming  flow  experiences  a  5% 
increase  in  static  pressure.  The  downstream 
boundary,  AY,  coincides  with  the  theoretical 

inviscid  shock  line.  Two  bleed  parameters 

(abbreviated  B.P.)  are  defined  in  order  to 
characterise  bleed: 


B.P  1 


(3) 


where  the  subscripts  w  and  oo  denote  plate 
and  free  stream  conditions  respectively,  u  la 
the  free  stream  velocity  of  the  incoming  flow 
and  the  v  denotes  the  vertical  velocity  at  the 
plata.  This  parameter  (B.P.  I)  directly 
specifies  the  mass  flux  at  each  grid  point  in 
tha  bleed  region.  The  flux  values  at  the 
boundaries  of  the  bleed  region  are  ramped  up 
to  the  full  value  over  typically  three  grid 
points  to  minimise  boundary  layer  effects  due 
to  sudden  changes  in  boundary  conditions.  It 
is  apparent  that  the  suction  under 
consideration  may  be  categorised  as  porous 
bleed  as  opposed  to  slot  bleed. 


Since  B.P  l  does  not  incorporate  the  area  of 
bleed  and  therefore  does  not  adequately 
describe  the  total  amount  being  bled  off,  a 
second  parameter  (B.P.  II)  is  defined: 
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where  I  =  Ag/l^,  M  is  the  total  mass  bleed 
rale  from  the  bleed  area  (A8)  and  BL  is  the 
mast  entering  the  domain  from  an  equal  area 
of  height  adjacent  to  the  plate  at  the 
upstream  boundary.  Four  cases  are  computed 


categorised  by  the  values  of  two  bleed 
parameters.  These  are  summarised  in  Table  2. 
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cm) 

p,JkP‘)  VO 

Expt. 

1.4 

2.93 

3.8x10s 

690  251 

Theory 

1.3 

2.93 

8.8x10s 

690  251 

Tabla 

X: 

Flow  Conditions  for  Experiment  and  Theory 

B.P.  I 

b.p  n 

Case  1 

0.000 

0.000 

Case  2 

0.010 

0.011 

Caa  3 

0.025 

0.029 

Case  4 

0.050 

0.055 

Table  3:  Bleed  Parameter  Values  Employed 


3.S  Details  of  Computations: 

Ths  computational  grid  employ*  32  streamwise  grid 
piano,  uniformly  spaced  in  the  x>direction  with  Ax  = 
tm.  The  upstream  boundary  is  located  at  54^  upstream 
of  the  fin  leading  edge,  and  the  downstream  boundary  at 
x  =  284^.  In  the  vertical  (y)  direction,  48  grid  points 
are  employed  integrating  up  to  approximately  84^  for 
Cases  1  and  2  and  8.54^  for  Cases  3  and  4.  The  number 
of  ordinary  point*  in  the  span  wise  direction  in  all  cases  is 
32.  The  sublayer  in  proximity  to  the  fin  is  resolved  with 
8  grid  points  for  all  cases  while  close  to  the  flat  plate  8 
pouits  are  employed  for  Cases  1  and  2  and  12  points  for 
Casa  3  and  4.  This  refinement  is  necessary  to  adequately 
resolve  all  relevant  flow  features.  A  complete  set  of 
criteria  for  judging  the  acceptability  of  a  grid  system  do 
not  exist  at  present.  A  number  of  necessary  criteria 
atablisbed  in  the  literature  are  presented  below. 

•  Viscous  Sublayer:  The  average  height  of  the 
first  grid  point  (located  within  the  sublayer) 
above  the  boundary  should  satisfy  [24|: 

*’*. 

g-*y  =  — ,v  <  j  (5) 
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V 

s.  -  -  0. 

and  N  is  the  bleed  correction  factor  of  Cebeci 
defined  earlier.  This  requirement  is  dictated  by 
the  use  of  an  algebraic  eddy  viscosity  model 
(cf.  16,  201.  In  the  presence  of  bleed,  the 
height  of  the  First  grid  point  must  also  satisfy: 


where  m  1  is  the  magnitude  of  bleed  imposed 
at  the  wall.  This  requirement  arisa  from  an 
analysis  of  the  sublayer  equations  in  the  region 
near  the  boundary  where  m  >  >  < 

•  Height  of  Sublayer  Region:  As  indicated 


> 


earlier,  the  average  height  of  the  sublayer 
region  is  restricted  to  17,  18i: 


Ti.  *  Tu 


where  the  jubecript  m  denotes  the  edge  of  the 
sublayer. 


Boundary  layer  A  minimum  of  fifteen  (15) 
grid  points  (including  sublayer  points)  are 
required  within  the  boundary  layer  [31,  32]. 


The  characteristics  of  the  grids  employed  for  each  of 
cases  satisfy  all  the  constraints  described  above  and  are 
presented  in  Table  3. 


Requirement 

r 

z-N 

*■ 

J 

NPBL 

Avg. 

Case  1 

1.47 

1.47 

47.41 

0.00 

16 

Max. 

3.03 

3.03 

97.41 

0.00 

Avg. 

Case  2 

1.47 

1.46 

50.4 

0.30 

16 

Max. 

3.01 

3.01 

97.01 

0.41 

Avg. 

Case  3 

0.96 

0.94 

43.1 

0.11 

21 

Max. 

3.07 

3.07 

96.8 

0.15 

Avg. 

Casa  4 

0.97 

0.93 

44.88 

0.23 

21 

Max. 

305 

3.05 

136.9 

0.3 

Legend:  NPBL  •  Number  of  Points  in  Boundary  Layer 


Tabla  3:  Grid  Characteristics 


Computations  are  carried  out  at  the  four-pipe  CYBER 
205  at  NASA  Ames  Research  Center.  A  typical  speed  of 
260-  MFLOPS  is  achieved  for  the  highly  vectorised 
portions  of  the  code.  The  flow  development  and  required 
CPU  times  are  provided  in  Table  4.  In  this  tabla,  Tc  is 
the  characteristic  time,  the  time  required  for  a  particle  at 
upstream  conditions  to  traverse  the  length  of  the 
computational  domain. 


Case 

T 

CPU  (hours) 

1 

3.5 

4.2 

2 

4.2 

3.6 

3 

5.1 

15.0 

4 

5.4 

17.5 

Tabla  4:  Flow  Development  and  CPU  hours 


2.4  Evaluation  of  Error  Associated  with  the 
Sublayer  Equations: 


The  effect  of  the  approximations  inherent  in  the 
sublayer  model  are  evaluated  by  integrating  the  full  3-D 
boundary  layer  equations  across  the  sublayer  region  (0  < 
>  <  tm  )  to  yield  the  following  general  expressions  for 
the  two  components  of  wail  shear  stress: 
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In  these  equations,  the  subscript  1  indicates  equivalence 
to  the  expression  employed  for  in  the  sublayer  model, 
the  subscript  2  indicates  the  correction  to  the  wall  shear 
stress  due  to  the  terms  neglected  in  the  sublayer 
equations  and  the  subscript  m  indicates  the  value  of  the 
quantity  at  the  edge  of  the  sublayer  The  axis  i’  is 
normal  to  the  surface  and  points  towards  the  interior  of 
the  domain,  the  axes  x’  and  y’  lie  in  the  plane  of  the 
surface  and  x’,  y’  and  s’  form  a  local  right-handed 
orthogonal  coordinate  system.  The  velocities  u',  v’  and  w' 
are  components  of  the  velocity  in  the  x’,  y'  and  s' 
directions  respectively.  The  quantity  m  is  the  bleed  mass 
flow  at  the  solid  boundary  and  is  assumed  to  be 
orthogonal  to  it.  The  error  terms  in  equations  (8)  and 
(9)  (TJm  and  TJy  respectively)  are  evaluated  with  the 
computed  solution.  It  may  be  noted  that  the  term  (;v'  - 
m)  is  assumed  identically  zero  at  all  grid  points  where 
bleed  is  applied.  An  error  measure  is  computed  as: 
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The  average  error  over  all  solid  boundaries  in  the 
domain  are  summarized  in  Table  5.  The  error  is 


comparable  to  the  accuracy  of  many  of  the  experimental 
measurements,  and  may  be  reduced  by  further  grid 
refinement.  It  may  be  noted  that  the  maximum  error  in 
skin  friction  associated  with  the  2-D  sublayer  equations 
applied  to  2-D  turbulent  interactions  was  observed  to  be 
approximately  0.2%  33 


Case  1 

2 

3 

4 

Average 

Error  (%)  5.1 

5.3 

2.4 

2.2 

Table  5:  Average  Error  Associated  With  Sublayer 

Approximation 


S  Results 

The  computed  results  are  evaluated  and  compared  with 
the  boundary  layer  measurements  of  Shapey  and 
Bogdonoff  10  and  the  surface  pressure  measurements  of 
Goodwin  6  .  It  is  emphasised  that  these  experiments  did 
not  impose  any  suction  and  are  presented  here  to 
highlight  the  effects  of  bleed.  In  the  following  discussion, 
the  upstream  boundary  layer  thickness  is  utilised  to 
scale  distances  (see  [9|  for  alternative  scaling  factors  for 
3-D  turbulent  interactions),  horisontal  distances  (xj  are 
measured  relative  to  the  location  of  the  inviacid  shock 
(x«kotk)  at  the  relevant  spanwise  position  and  flow 
variables  are  normalized  by  the  value  of  the  variable 
under  free  stream  conditions  unless  otherwise  mentioned. 
The  experimental  stations  of  Shapey  are  located  as  shown 
m  rig.  3.  The  coordinates  of  each  location  are  provided 
id  Table  6. 


Survey  Location 

xJ*m 

*/*. 

1 

-5.40 

5.81 

2 

-4.40 

5.81 

3 

-3.40 

5.81 

4 

■2.40 

5.81 

5 

-1.40 

5.81 

6 

-0. 40 

5.81 

7 

0.60 

5.81 

8 

2.60 

5.81 

9 

3.94 

7.81 

10 

-0.14 

4.81 

11 

1.13 

3.81 

Legend:  x,=x  -  x,^ 


Table  6:  Survey  Locations  for  the  Experiments  of 

Shapey  (no  bleed) 

3.1  Effect  of  Suction  on  Flow  Variables: 

The  effect  of  suction  on  pitot  pressure  profiles  at 
Stations  1  and  2,  located  upstream  of  the  interaction  and 
at  the  experimentally  determined  line  of  upstream 
influence  respectively  is  negligible.  These  profiles  are 
similar  to  those  observed  in  2-D  boundary  layer  flows. 
At  Station  3  (not  shown),  located  in  the  vicinity  of  the 
line  of  coalescence,  the  pitot  pressure  profiles  for  Cases  2 
to  4  are  slightly  lower  than  Case  1  up  to  about  y  = 

1  beyond  which  the  profiles  .match  precisely.  At  Stations 
4  (not  shown).  5  (Fig.  4)  and  6  (not  shown),  located 
between  the  line  of  coalescence  and  the  theoretical 
inviscid  shock  location,  the  previously  observed  10 
overshoot  is  significantly  reduced.  An  approximately  30% 
reduction  in  pitot  pressure  is  observed  at  y  6  '  1 
between  Cases  1  and  4  at  Station  5  This  overshoot  is 
associated  with  the  compression  system  ahead  of  the 
shock  wave,  the  effect  of  suction  on  which  is  discussed 
later  At  Stations  7  (not  shown)  and  8  (Fig.  5)  located 
immediately  downstream  of  the  shock,  the  numerical 


model  does  not  display  the  prominent  overshoot  observed 
experimentally.  The  slight  S-shaped  behavior  near  the 
plate  is  more  apparent  at  these  stations.  It  is  clear  from 
the  pitot  pressure  profiles  presented  that  the  effect  of 
bleed  is  to  increase  the  pitot  pressure  modestly  up  to 
distances  of  y/<M  approximately  equal  to  1.  At  roughly 
this  distance  from  the  flat  plate,  the  curves  cross 
indicating  lower  pitot  pressure  values  in  the  presence  of 
suction  for  y ft  greater  than  approximately  1.  The 
pitot  pressure  profile  at  Station  10  (not  shown)  shows 
behavior  qualitatively  similar  to  that  observed  at  Station 
7.  The  lowest  pitot  pressure  value  at  Station  10  is 
overpredicted  numerically  by  about  25%  (for  the  no  bleed 
case)  at  y /<  *  0.4.  Station  9  displays  behavior  similar 
to  Station  3  and  is  not  shown. 

The  effect  of  suction  on  the  yaw  angles  is  now 
compared  with  the  experimental  data  of  Shapey.  At 
Stations  1  and  2  (not  shown),  no  effect  of  suction  is 
observed  on  the  yaw  angles  which  are  relatively  small.  At 
Station  3  (Fig.  6),  significant  reduction  is  observed  in  the 
yaw  angle  in  the  region  up  to  y/f  <  0.8.  The  small 
negative  yaw  angle  observed  for  Case  4  is  currently  under 
investigation  and  is  thought  to  be  due  to  the  close 
proximity  of  this  station  to  a  grid  point  where  boundary 
conditions  dictate  significant  bleed.  Stations  4  (not 
shown)  and  5  (Fig.  7)  also  display  modestly  lower  vaw 
angles  for  larger  values  of  bleed.  At  Station  6  (Fig.  8), 
a  'crossover’’  of  the  profiles  is  observed  similar  to  that 
observed  for  pitot  pressure.  Yaw  angles  are  smaller  with 
higher  suction  near  the  fiat  plate  and  larger  at  distances 
further  away.  Similar  results  are  observed  at  Stations  7 
through  11  (not  shown). 

The  effect  of  bleed  on  the  compression  system  upstream 
of  the  shock  is  more  clearly  visible  in  contour  plots  of 
static  pressure.  Figs.  9  through  12  show  pressure 

contours  at  a  streamwise  location  of  1 14  downstream  of 

30 

the  fin  leading  edge  for  the  four  cases  respectively  The 
shock  wave  is  located  in  the  region  of  most  contour 
density  and  is  clearly  visible  in  the  inviscid  region  A 

perceptible  increase  in  the  concentration  of  pressure 
contours  forming  the  nght  ’leg”  of  the  shock  wave  is 

evident  leading  to  the  conclusion  that  suction  has  the 

effect  of  tightening  the  compression  fan  It  may  be  noted 
that  these  figures  also  provide  an  indication  of  the  extent 
of  shock  smearing  due  to  the  numerical  algorithm. 

The  computed  wall  pressure  is  compareq  with  the 

experimental  data  of  Goodwin  61  in  Fig  13  in  which  the 
surface  pressure  is  plotted  against  streamwisp  distance 
(measured  relative  to  the  theoretical  inviscid  shock 
location  xth  at  a  spanwise  location  of  r  ^  =  6  9 

This  figure  indicates  that,  for  the  no  bleed  case,  the 
computations  accurately  predict  the  extent  of  upstream 
influence  and  the  pressure  rise  in  the  interaction  region. 
The  effect  of  bleed  is  evidently  to  retard  the  point  of 

upstream  influence  slightly  in  the  downstream  direction 
thus  reducing  moderately  the  total  distance  over  which 
the  pressure  rise  is  achieved  It  may  be  noted  that  the 

experimental  data  of  Hingst  and  Tanji  34  for  2-D 
turbulent  Interactions  in  the  presence  of  bleed  display 
qualitatively  similar  effects.  The  overall  pressure  rise  :s 
not  affected  in  any  significant  manner  The  pronounced 
drop  in  wall  pressure  (Fig  13)  at  x-x  approximates 
equal  to  -44^  is  probably  due  to  relatively  large  iocai 
values  of  bleed  that  exist  in  this  region  for  Cases  3  and 
4 


The  localised  effect  of  bleed  i a  indicated  in  the  surface 
skin  friction,  shown  in  Fig.  14  at  the  same  location  aa 
the  data  of  Goodwin.  A  steep  rise  in  skin  friction  is 
observed  in  the  bleed  region  with  negligible  differences 
elsewhere.  The  peak  values  Increase  with  suction.  It  is 
apparent  that  there  is  a  significant  drag  penalty 
associated  with  the  employment  of  suction. 

Eddy  viscosity  values  are  observed  to  be  generally  lower 
in  the  presence  of  bleed.  The  eddy  viscosity  at  station  5 
(normalised  by  the  free  stream  dynamic  viscosity,  )ij  is 
plotted  in  Fig.  15  at  Station  5  and  shows  an 
approximately  70%  drop  in  value  between  Cases  1  and  4 
at  y/j"  =  0.4.  Previous  research  has  indicated  however, 
that  the  details  (e.g.,  velocities,  pressures  and 
temperatures)  of  3-D  turbulent  interactions  are  relatively 
insensitive  to  the  particular  turbulence  model  employed 
with  the  exception  of  a  small  fraction  of  the  boundary 
layer  adjacent  to  the  wall. 

3.3  Effect  of  Suction  on  Flow  Field  Structure: 

The  computed  solutions  are  utilised  to  examine  the 
effect  of  bleed  on  the  flow  structure.  Figs.  16  through  19 
show  computed  surface  skin  friction  lines  (*  surface 
streamlines*)  for  Cases  1  through  4  respectively  The 
location  of  the  inviacid  shock  wave  is  indicated.  It  is 
evident  that  the  lines  of  coalescence  and  divergence 
persist  in  all  cases.  The  line  of  coalescence  indicates  a 
tendency  to  align  with  the  shock  location  in  the  limiting 
case,  itus  limiting  position  is  probably  dictated  by  the 
particular  region  of  bleed  under  consideration  which 
terminates  at  the  invisdd  shock  line. 

Further  investigation  of  the  computed  flow  fields  is 
carried  out  with  particle  tracing  methods.  A  number  of 
particles  are  released  in  a  systematic  fashion  at  the 
upstream  leading  edge  of  the  domain  and  their  motion  is 
numerically  integrated.  Particles  in  sets  of  ten  released  at 
two  heights  above  the  flat  plate  are  presented  to  bring 
out  the  effects  of  suction  on  flow  structure.  Two  views 
are  shown  for  each  case  and  each  height  -  a  top  view 
and  a  view  looking  downstream.  The  flat  plate  is  outlined 
in  each  figure  for  clarity. 

Figs.  20  through  23  show  traces  of  particles  released  at 
a  height  of  y/f —  =  0.1.  Particles  for  Casas  1  and  2  enter 
the  vortical  structure  with  no  obvious  differences.  For 
Case  3  on  the  other  hand,  whereas  particles  closer  to  the 
fin  enter  the  vortex,  those  away  from  the  fin  are  in  fact 
ingested  into  the  flat  plate.  This  ingestion  is  even  more 
apparent  for  Case  4  where  the  bleed  is  highest  and  all 
but  one  particle  -  that  closest  to  the  fin  -  are  ingested. 
It  may  be  mentioned  that  near  the  corner  formed  by  the 
fin  leading  edge  and  the  flat  plate,  the  particle  traces 
may  be  inaccurate  due  to  the  large  gradients  present.  It 
is  clear  that  with  increase  in  suction  a  larger  portion  of 
the  boundary  layer  is  ingested  as  may  be  expected. 

Higher  into  the  boundary  layer,  all  cases  display  the 
prominent  vortical  structure  mentioned  previously  Figs. 
24  through  27  show  traces  of  particles  released  at  y  ■  j 
=  0.5.  For  ail  cases,  particles  nearer  the  fin  are  swept 
underneath  those  further  away.  The  particles  rotate 
counterclockwise  when  viewed  in  the  downstream 

direction.  It  is  apparent  that  the  effect  of  bleed  on  the 
shape  and  dimension  of  the  vortical  structure  is 

negligible.  This  is  more  obvious  in  Fig.  28  where  the 


trace  of  a  single  particle  (originating  at  a  spanwise 
distance  of  5 and  a  vertical  distance  of  0.54^  above 

the  flat  plate)  is  shown  for  all  cases.  At  this  distance 

above  the  flat  plate  (y/4^)  =  0.5,  the  yaw  angle  profiles 
discussed  previously  (see  Figs.  6  through  8)  display 

modestly  lower  values  In  the  presence  of  bleed.  The 
particle  traces  in  Fig.  28  are  therefore  consistent  with  the 
yaw  angles  presented  previously. 

It  is  clear  that  the  mean  flow  field  pattern  is 

dominated  by  the  large  vortical  structure  even  in  the 
presence  of  strong  bleed.  The  effect  of  bleed  is  to  simply 
ingest  a  larger  portion  of  the  boundary  layer  for  larger 
suction  values. 


4  Weak  Shock-Shock  Intersection  in  the 
Presence  of  a  Turbulent  Boundary  Layer 

Mee  [35]  and  Mee  and  Stalker  [15]  describe  several 
experiments  with  weak  shock-shock  intersection  in  the 
presence  of  a  turbulent  boundary  layer.  They  conclude 
that  ’intersecting  shock  interactions  can  produce  a  given 
overall  pressure  rise  with  less  likelihood  of  separation 
than  an  equivalent  strength  single  shock  interaction*.  As 
a  benchmark,  the  flow  past  two  5*  fins  mounted  on  a 
flat  plate  (Fig.  29)  is  computed  and  compared  with 
experiment.  The  governing  equations  and  numerical  model 
are  as  described  in  Section  2.1.  The  flow  parameters  are: 

*  1  85 
»  1.9  (mm) 

/fc-4.1  107  (m“ l) 

Pt  -  S00  [kPa) 

T~  -  295  (A) 

*00 

4>t  is  the  average  boundary  layer  thickness  in  the 
experimental  measurement  region  jlSj.  Since  line  EF 
(Fig.  29)  is  a  line  of  symmetry,  the  flow  in  the  half  area 
ABFE  is  computed.  The  flow  in  EFCD  is  then 
constructed  by  symmetry  considerations.  On  the  fin  and 
plate  surfaces,  solid  boundary  conditions  are  applied  (see 
fin  surface  in  Section  2.2).  At  the  upstream  boundary, 
the  flow  is  assumed  to  be  2-D  and  as  for  the  bleed  cases, 
the  flow  variables  are  generated  with  a  separate  2-D  code 
by  matching  the  average  boundary  layer  thickness  with 
that  published  in  the  experiments  of  Mee  and  Stalker. 
Zero  gradient  extrapolation  is  emploved  at  the 
downstream  boundary  and  symmetry  boundary  conditions 
are  applied  on  the  symmetry  boundaries.  The  grid 
employs  forty  (40)  streamwise  grid  planes,  uniformly 
spaced  in  the  x-direction  with  ,ix=4tT.  The  upstream 
boundary  is  located  44>t  upstream  of  the  fin  leading  edge, 
and  the  downstream  boundary  is  at  354t>  downstream. 
Forty  eight  (48)  and  thirty  two  (32)  grid  points  are 
respectively  employed  in  the  vertical  (y)  and  spanwise  (r) 
directions.  Both  the  fin  and  plate  sublayers  are  resolved 
with  eight  (8)  points  each.  All  grid  constraints  described 
in  Section  2.3  are  satisfied  and  details  are  omitted  for 
the  sake  of  brevity  Convergence  is  achieved  at  a  flow 
development  time  of  5.2  T  (see  Section  2.3) 
corresponding  to  a  CPL'  requirement  of  5.5  hours  on  the 
CYBER  205 


The  available  experimental  results  35  consist  of  surface 
pressure  measurements  in  the  interaction  region  and 
limiting  streamline  patterns  recorded  with  the  china  film 
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technique.  Fig.  30  shows  numeric*]  pressure  contours  on 
the  flat  plate.  The  quantity  plotted  in  this  figure  is  the 
pressure  rise  above  its  upstream  value  normal i ted  with 
the  theoretical  inviscid  pressure  rise  due  to  the  first  shock 
wave.  The  experimental  surface  pressure  contour  plot  is 
shown  in  Fig.  31  and  corresponds  to  the  area  marked 
AHXY  in  Fig.  30.  A  comparison  of  the  two  figures 
reveals  differences  in  the  extent  of  upstream  influence  of 
the  interaction  and  the  rate  of  pressure  rise.  The 
computations  predict  a  larger  upstream  influence  than 
observed  in  experiment.  This  is  in  contrast  to  results 
published  previously  [10]  where,  for  stronger  (single  shock) 
interactions  at  higher  Mach  and  Reynolds  numbers, 
computations  with  two  different  eddy  viscosity  models 
both  underpredicted  the  extent  of  upstream  influence.  In 
the  present  instance,  if  the  upstream  influence  is 
measured  as  the  distance  normal  to  the  shock  wave  where 
the  pressure  rises  by  10%  of  the  final  pressure  rise  after 
the  first  shock,  the  computed  upstream  influence  distance 
is  approximately  13%  larger.  The  computed  surface 
pressure  also  increases  at  a  steeper  rate  after  the  initial 
pressure  rise.  Fig.  32  shows  computed  surface  streamlines 
and  may  be  compared  with  those  obtained  experimentally 
(Fig.  33).  The  resemblance  of  the  two  patterns  - 
numerical  and  experimental  •  to  each  other  is  clearly 
evident.  As  expected  for  such  weak  interactions,  there  is 
no  line  of  coalescence.  Further  analysis  of  the  computed 
flow  field  indicates  the  absence  of  separation  (as 
expected).  Fig.  34  shows  traces  of  particle*  released  at 
0.1  and  0.3  t n  above  the  flat  plate  respectively. 

5  Conclusion*  and  Future  Work 

A  computational  study  of  the  effect  of  bleed  on  the  3-D 
shock  wave  turbulent  boundary  layer  interaction  generated 
by  a  sharp  fin  at  Mach  3  for  fin  angle  a  =  20°  and 
Rej  =  9  x  10*  indicates  remarkably  modest  influence  of 

suction.  The  effect  of  bleed  is  to  cause  a  tightening  of 
the  compression  system  ahead  of  the  shock  wave.  The 
overshoot  observed  in  pitot  pressure  profiles  due  to  the 
compression  system  ahead  of  the  shock  is  significantly 
reduced.  There  is  a  tendency  of  the  line  of  coalescence  to 
align  with  the  downstream  boundary  of  bleed  although 
this  statement  is  probably  valid  only  for  the  bleed  region 
under  present  consideration.  The  computations  indicate 
that  the  particular  bleed  schedule  described  is  incapable 
of  controlling  or  significantly  modifying  the  undesirable 
separation  and  large  vortical  structure  observed  in 
previous  research. 

Results  from  a  numerical  simulation  of  a  weak  shock- 
shock  intersection  in  the  presence  of  a  turbulent  boundary 
layer  show  overall  good  agreement  with  available 
experimental  data.  Current  research  efforts  are  focussed 
on: 

1.  Bleed  behind  shock:  For  completeness,  the 
effect  of  1%  bleed  applied  in  the  area  between 
the  fin  and  the  theoretical  inviscid  shock  wave 
is  being  investigated. 

2.  Shock-Shock  Interactions:  It  is  proposed  to 
examine  the  flow  structure  due  to  a  strong 
shock-shock  intersection  in  the  presence  of  a 
turbulent  boundary  layer.  The  precise 
parameters  to  be  employed  are  currently  being 
determined  in  conjunction  with  the  Princeton 
Gas  Dynamics  Laboratory 
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